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ABSTRACT 

An experimental investigation was conducted to determine the effect of various 
physical parameters on the pressurant gas requirements during the pressurization and 
expulsion of liquid hydrogen from a 1. 52-m- (5 -ft-) diameter spherical tank. The exper- 
imental results were compared with results predicted by an analytical program. Tests 
were conducted for a range of outflow rates at a nominal operating pressure level of 
34. 47xl0 4 N/m 2 (50 psia) using inlet gas temperatures of 167, 298, and 355 K (301°, 
536°, and 639° R). Data were obtained using three pressurant injector geometries 
(hemisphere, radial, and straight pipe) and two tank wall thicknesses (0. 762 and 
0. 409 cm; 0. 30 and 0. 161 in. ). The experimental results indicate that the inlet gas 
temperature has the strongest influence on the pressurant gas requirements followed 
closely by the influence of injector geometry. Fair agreement existed between the anal- 
ysis and experimental data when using a diffuser -type (hemisphere and radial) injector. 
The agreement between the experimental data and the analysis was not good when using 
the straight pipe injector. 
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SUMMARY 

An experimental investigation was conducted to determine the effect of various 
physical parameters on the pressurant gas requirements during the pressurization and 
expulsion of liquid hydrogen from a 1. 52-meter- (5 -ft-) diameter spherical tank. Tests 
were conducted for a range of outflow rates at a nominal operating pressure level of 
34. 47x10^ newtons per square meter (50 psia) using inlet gas temperatures of 167, 298, 
and 355 K (301°, 536°, and 639° R). Data were obtained using three pressurant injector 
geometries (hemisphere, radial, and straight pipe) and two tank wall thicknesses (0. 762 
and 0. 409 cm; 0. 30 and 0. 161 in. ). 

The experimental results indicate that the inlet gas temperature has the strongest 
influence on the pressurant gas requirements followed closely by the influence of injector 
geometry. Results also indicate that tank wall thickness does not greatly influence the 
pressurant gas requirements. 

The experimental results were compared with results predicted by a previously 
developed analytical program which was revised and extended for the test conditions. 

Fair agreement existed between the analysis and experimental data when using a diffuser- 
type injector. The agreement between experimental data and the analysis was not as good 
when using the straight pipe injector. 


INTRODUCTION 

During the past several years, a great deal of effort has been devoted to the prob- 
lems associated with the pressurized discharge of a cryogenic liquid from a tank. The 
main objectives of these efforts have been toward the optimization of a propellant tank 
pressurization system. One phase of this optimization is a precise determination of 
pressurant requirements for any given set of operating parameters (e.g. , tank pressure, 


inlet gas temperature, liquid flow rate, etc.)- This knowledge would allow the design of 
a system that carried only the weight of gas necessary to accomplish the mission. 

Several investigators have developed analyses (e.g. , refs. 1 and 2) which attempt to 
predict the pressurant gas requirements during the pressurized discharge of a cryogenic 
fluid. These analyses, however, are either burdened with simplifying assumptions or 
involve parameters and terms about which little is generally known "a priori. ” Because 
of these limitations the validity of the analytical results has to be verified largely by 
correlations of experimental results. 

Therefore, an investigation was conducted at the Lewis Research Center to deter- 
mine experimentally the effect of various physical parameters (e.g. , tank geometry and 
size, injector geometry, pressurant gas, inlet gas temperature, etc. ) on the pressurant 
gas requirements during the expulsion of liquid hydrogen from oropellant tanks. The 
information obtained could be used in the design of optimum pressurization systems for 
particular missions. The results of this investigation were also used to determine the 
capability of the analysis of reference 1 to predict the pressurant gas requirements for 
the various physical parameters. Although the analysis of reference 2 is widely used in 
the design and analysis of pressurization systems, its use requires values for many 
empirical constants and physical parameters. In many cases these values must be 
obtained from experimental data on prototype systems. Inasmuch as these constants 
and parameters were not available, it was decided to use the analysis of reference 1 
since it required the least amount of prior knowledge of a particular system. 

Small scale test work was done in a 0. 82 cubic meter (29 ft ) cylindrical tank to 
determine the effect of injector geometry on the pressurant gas requirements (GHg) 
during the dischange of liquid hydrogen (LH 2 ). The results of this work are reported in 
reference 3. The analysis, for this case, proved to be adequate (within ±10 percent) in 
predicting the pressurant gas requirements even though two of its major assumptions 
were shown to be invalid - namely, no heat transfer to the liquid surface, and no mass 
transfer . 

This report presents the results obtained in 1. 52 -meter- (5-ft-) diameter spherical 
aluminum tanks. The primary objective of these tests was to obtain additional experi- 
mental information on the pressurant requirements during the initial pressurization as 
well as the expulsion period in larger scale tanks for the various operation parameters. 
Information was also obtained on tank wall heating, liquid heating, residual ullage energy, 
and mass transfer in order to gain an insight into the reasons for any variations in pres- 
surant gas requirements. An additional objective was to compare the experimental 
results with the predicted ones. 

For this report the analysis of reference 1 was revised and extended (as described 
in appendixes A and B) to account for energy transfer occurring at the gas -liquid interface 
in tanks of any arbitrary symmetric tank shape. The analysis was also modified and 
extended (as described in appendix C) to cover the initial pressurization period. The 
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major limiting assumptions still remaining from the original analysis are one- 
dimensional flow and no mass transfer. Although reference 3 showed tie inadequacy 
of these assumptions when analyzing a straight pipe injector, a more ri orous analysis 
was not attempted because of the added complexity. 

The tests reported herein were conducted in a vacuum chamber using 1. 52 -meter - 
(5— ft-) diameter spherical aluminum tanks. The main test variables were as follows: 

Tank wall thickness, cm (in. ) 0. 762 and 0. 409 (0. 30 and 0. 161) 

Pressurant gas injector geometry hemisphere, radial, and straight pipe 

Nominal inlet gas temperature, K (°R) 167, 298, and 355 (301, 536, and 639) 

Liquid outflow rates, kg/ sec (lb/ sec) between 0. 859 and 0. 218 (1. 890 and 0. 481) 

Initial gas ullages in tank, percent 4, 28, 55, and 75 

The high outflow rate corresponds to a tank expulsion time of 132. 4 seconds and a thrust 
level of 2290 kilograms (5050 lb). The expulsion time and thrust level corresponding to 
the low flow rate are 531. 9 seconds and 580 kilograms (1280 lb). The thrust levels are 
based on a specific impulse of 444 seconds and an oxidant to fuel ratio (O/F) of 5. Data 
were also obtained for pressurization of the tank at various rates (6. 48x10^ to 2. 65x10^ 

o 

(N/m )/sec; 0. 94 to 3. 85 psi/sec) from 1 atmosphere to the desired operating level. All 
tests were performed using gaseous hydrogen as a pressurant at a nominal operating tank 
pressure of 34. 47x10^ newtons per square meter (50 psia). 


SYMBOLS 


A 

b 

C 

C 


H 


w 


area, m^ (ft^) 


1 - 


2 — 


Z r 

orifice coefficient 

effective perimeter of interior hardware, m (ft) 
specific heat at constant pressure, J/ (kg)(K) (Btu/ (lb)(°R)) 
specific heat at constant volume, J/(kg)(K) (Btu/(lb)(°R)) 
specific heat of tank wall, J/ (kg)(K) (Btu/ (lb)(°R)) 


a - (o!Ci>) - T. 


w 


\ w^w/ 


| (At) — 


'w 


D 


orifice diameter, m (ft) 


3 



p p 

g gravity acceleration, m/sec (ft/ sec ) 

H enthalpy, J (Btu) 

h specific enthalpy, J/kg (Btu/lb) 

h convective heat -transfer coefficient, j/(m 2 )(K)(sec) (Btu/(ft 2 )(°R)(sec)) 

k thermal conductivity, j/(m)(K)(sec) (Btu/(ft)(°R)(sec)) 

L flow length, m (ft) 

l thickness, m (ft) 

M mass, kg (lb) 

AM differential mass, kg (lb) 

M mass flow rate, kg/ sec (lb/sec) 

M molecular weight, kg/ (kg) (mole) (lb/ (lb) (mole)) 

Mqn net mass added to ullage, Mq± M t , kg (lb) 

Mj ideal pressurant requirement, kg (lb) 

N number of volume segments 

Nu Nusselt number, h c L/k 

n (or i) summing index 
P pressure, N/m 2 (lb/in. 2 ) 

Pr Prandtl number, pCp/k 

p p 

Ap differential pressure, N/m (lb/in. ) 

Ap* orifice Ap 

Q heat transfer, J (Btu) 

Q heat -transfer rate, J/sec (Btu/sec) 

Q’ specific heat-transfer rate, J/ (kg)(sec) (Btu/ (lb)(sec)) 

q heat-transfer rate per unit area, J/ (m )(sec) (Btu/ (ft )(sec)) 

R gas constant, J/ (K)(mole) 

Re Reynolds number, pLV//i 
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r 

Ar 

T 

AT 

T 6 

t 

At 

U 

AU 

u 

V 

AV 

V 
v 
W 

*n 

x 

Ax 

V 

y 

z 


a 


0 

v 

6 


radius, m (ft) 
increment of radius, m (ft) 
temperature, K (°R) 
differential temperature, K (°R) 

temperature at the edge of thermal boundary layer, K (°R) 
time, sec 

time increment, sec 

internal energy, J (Btu) 

differential energy, J (Btu) 

specific internal energy, J/kg (Btu/lb) 

volume, m** (ft**) 

volume increment, m** (ft**) 

velocity, m/sec (ft/sec) 

specific volume, m /kg (ft /lb) 

work, J (Btu) 

number of volume elements in the ullage 
coordinate in direction of tank axis, m (ft) 
space increment, m (ft) 
expansion factor 

thickness within the boundary layer, m (ft) 

compressibility factor, Z = pv/RT 

elevation or vertical distance along tank wall, m (ft) 


1 + 


h c At 

U w P w ^wy 


(2h c RZ At> 
rMCP 

r 


1 + 



1/2 


coefficient of thermal expansion, 1/K (1/°R) 
specific heat ratio 

finite increment, or total boundary layer thickness, m (ft) 
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■Ill 


X 

M 

P 

o> 


latent heat of vaporization, J/kg 
viscosity, kg/ (m)(hr) (lb/(ft)(hr)) 
density, kg/m^ (lb/ft 

(— Z 1 + 

\ M il Mjrr 



(Btu/lb) 


Subscripts: 


A analytical results 

ad adiabatic 

E experimental results 

f final state or condition 

G gas added to tank 

H internal hardware 

1U 

i initial state or condition, or i m term 

i— f from initial to final state or condition 

L liquid 

o condition prior to ramp 

S liquid surface 

sat saturation 

T total quantity 

t transferred 

U ullage 

w wall 


Superscript: 

1 time index 


APPARATUS AND INSTRUMENTATION 
Facility 

All tests were conducted under vacuum inside a 7. 61-meter- (25-ft-) diameter spher- 
ical vacuum chamber (fig. 1) to reduce the external convective heat transfer into the pro- 
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Figure 1. - 7. 61 -Meter- (25-ft-) diameter vacuum chamber. 
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pellant tank. The vacuum capability of this chamber was approximately 8x10 millimeter 
of mercury. A general schematic of the test tank and associated equipment is shown in 
figure 2. A heat exchanger and blend valve subsystem capable of delivering gaseous 
hydrogen at temperatures of 167 to 405 K (301° to 729° R) at a maximum flow rate of 

o 

5. 40xl0 - ' 6 kilogram per second (0. 12 lb/sec) were used to control pressurant gas inlet 
temperature. A ramp generator and control valve were used for controlling the initial 
rate of pressurization of the propellant tank. A closed-loop pressure control circuit was 
used to maintain constant tank pressure during the expulsion period. The liquid outflow 
rate was controlled by remotely operated variable flow valves. The liquid hydrogen out- 
flow from the tank was returned to the storage Dewar. 

Liquid outflow rates were measured using a turbine -type flowmeter located in the 
transfer line. The flowmeter was calibrated (within an estimated uncertainty of ±1/2 per- 
cent) with liquid hydrogen over the expected range of flow rates. The calibration was per- 
formed at Lewis Research Center. Pressurant gas inlet flow rates were determined by 
the use of an orifice located in the pressurant supply line. Tank, line, and differential 
pressures were measured with bonded strain -gage -type transducers (estimated uncer- 
tainty of ±1/4 percent). 
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Gaseous 



valves 


Figure 2. - General schematic of facility. 


Test Tanks 


The experimental work was conducted in 1. 52-meter- (5-ft-) diameter bare wall 
spherical al um i num tanks. One tank (fig. 3) had an average wall thickness of 0. 762 centi- 
meter (0. 30 in. ). The other tank was identical except the wall was chem-milled down to 
an average thickness of 0. 409 centimeter (0. 161 in. ). One stainless steel lid served both 
tanks. The lid housed the inlet and vent pipes and the electrical connections for all inter- 
nal tank instrumentation. The lid was 0. 457 meter (18 in. ) in diameter, 3. 18 centimeters 
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Figure 3. - 1.52-Meter- (5-ft-) diameter spherical test tank. Tank wall thickness, 
0.762 centimeter (0.30 in.). 


(1. 25 in. ) thick, and weighed 67. 50 kilograms (149 lb). The inner surface of the lid 
conformed to the contour of the tank wall. 

A view port and television camera were installed on the 0. 7 62 -centimeter (0. 30-in. ) 
wall test tank (fig. 4) for the purpose of visually locating the initial liquid in the tank prior 
to expulsion and to enable observation of any physical processes occurring in the tank. 
Lighting the tank interior was accomplished by bulbs mounted on the inner surface of the 
tank wall. 


Pressurant Gas Injector Geometries 


The three injector geometries tested are shown in figure 5. These particular geom- 
etries were selected because they represent the three basic injection patterns of refer- 
ence 3. The hemisphere injector diffuses the pressurant uniformly in all directions into 
the ullage volume. The radial injector injects the pressurant gas radially toward the 
tank wall. And the straight pipe injects the oressurant gas in a concentrated jet toward 
the liquid surface. The hemisphere and radial injectors, which best seem to conform to 
the analytical assumption of one -dimensional gas flow as shown in reference 3, have 
approximately the same open exit area for the hemisphere (176. 8 cm ; 27. 4 in. ) and for 
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Figure 4. - View port and television camera installation on test tank. Wall 
thickness, 0.762 centimeter (0.30 in.). 


— - Denotes pressurant gas flow paths 



(a) Hemisphere injector. Open 
area, 176.8 square centimeters 
(27. 4 in. 2). 


(b) Radial injector. Open area, 
170.3 square centimeters 
(26. 4 in. 2). 


(c) Straight pipe injector. 
Open area, 9. 5 square 
centimeters (1.47 in. 2 ). 


Figured - Injector geometries. All dimensions are in centimeters (in.). 
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2 2 

the radial injector (170. 3 cm ; 26. 4 in. ). The straight pipe injector had an open area of 

2 

9.5 square centimeters (1. 47 in. ). 

Internal Tank Instrumentation 

Ullage gas temperatures were used to determine the mass and energy content of the 
tank ullage. These temperatures must be obtained with sensors capable of accurate 
measurement of rapid changes in temperature. Internal tank instrumentation is illustrated 
in figure 6. Location of the vertical and horizontal ullage gas temperature rakes are 
indicated. The thermopile was the basic temperature measurement technique used in 
this investigation. The use of thermopiles to measure ullage gas temperature was first 
developed in reference 4, and the technique was used with good results in reference 3. 

The main advantage of using thermopiles is their fast response time (between 0. 2 and 1. 0 
sec) in going from saturated liquid to saturated vapor. This time response is approxi- 
mately an order of magnitude less than sensors generally used in this type of investi- 
gation. 

A typical three -element thermopile unit and its associated wiring schematic are 
illustrated in figure 7(a). The thermopile units were constructed of 0. 202 millimeter 
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Constantan 

Chrome! 

Copper 

Thermocouple 
junction (typical) 



Measurement level 
Reference level 


(a) Three element thermopile. 


fTypical 



resistor probe 

(b) Wiring schematic for rake. 



Figure 7. -Thermopile rake. 


chromel and constantan wire. The support structure was made of Bakelite to minimize 
thermal conduction from one measurement station to another. Vertical ullage gas temper- 
ature profiles were obtained by stacking the individual thermopile units as shown in fig- 
ures 7(b) and (c). The spacing between the reference and measuring levels was 2. 54 cen- 
timeters (1. 0 in. ) for the top 45 thermopiles composing the vertical rake. The 10 units 
at the bottom of the rake had a 1. 27-centimeter (0. 50-in. ) spacing in order to obtain a 
more accurate temperature profile of the ullage gas near the liquid surface at the end of 
an expulsion. 

Platinum resistance sensors, which were located at least every tenth station starting 
from the bottom of the rake, sensed the absolute temperature at their location and pro- 
vided a reference for the thermopiles above the location. The horizontal rakes were 
composed of platinum resistance sensors spaced a maximum of 12.70 centimeters (5. 00 
in. ) apart in the radial direction. Two platinum resistance sensors were used at each 
location to measure liquid and/or gas temperatures in the ranges 20 to 38. 9 K (36° to 
70° R) and 38. 9 to 277. 8 K (70° to 500° R). These dual sensors permitted more accurate 
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measurements of liquid and gas temperatures than could be achieved with one sensor 
covering the entire range. 

The initial static temperature profile near the liquid surface was determined by a 
fixed interface rake. This rake contained 13 platinum resistance sensors spaced 0. 64 
centimeter (0. 25 in. ) apart. The range of these sensors was 20 to 38. 9 K (36° to 70° R). 
For complete tank expulsions (i. e. , 5 percent ullage to 5 percent outage), the initial 
liquid level was always assured to be on the interface rake in the thick walled tank by 
monitoring a closed-circuit television system which viewed, through a port in the tank, 
the position of the liquid-gas interface on a scale fixed to the interface rake (fig. 6). For 
tests conducted at the 28, 55, and 75 percent ullage values, the scale was moved lower in 
the tank and the television camera was refocused. The liquid level for the thin-walled 
tank was determined by hot-wire liquid level sensors located on the interface rake. 

Platinum resistance sensors were also used to determine tank wall temperatures at 
twelve locations and the liquid temperature at the flowmeter. Copper -constantan thermo- 
couples were used to determine tank lid temperatures at five locations and the pressurant 
gas inlet temperature. 

All measurements were recorded on a high-speed digital data system. The measure- 
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ments were recorded at a rate of 3. 125x10 channels per second. Each measurement 
channel was sampled every 0. 064 second. 


PROCEDURE 

The spherical test tank was filled from the bottom to approximately a 2 percent ullage 
condition. It was then topped off as necessary while the tank lid and peripheral support 
hardware reached steady -state operating temperatures. 

Temperature conditioning of the pressurant gas was then started. Gas flow was 
established through the heat exchanger loop, through the control valves and orifice 
arrangement, and then into the tank ullage from which it was vented through the condi- 
tioning line to the outside as shown in figure 2 . The temperature control circuit shown in 
figure 2 was used to get the desired pressurant gas temperatures during the flow period. 
When the pressurant gas temperature conditioning was almost completed, the liquid level 
in the test tank was adjusted to a desired value by either topping or slow draining. The 
liquid level was observed on closed -circuit television for the thick -wall tank. Liquid 
level, for the thin-wall tank, was located by the use of hot-wire liquid level sensors. The 
pressurant gas flow was then stopped and the test tank was vented in preparation for an 
expulsion run. The automatic controllers and timers were preset with all the desired run 
and operating conditions (i.e. , tank pressure level, length of ramp period, length of hold 
period, liquid outflow valve position, start and stop times of the data recording equipment, 
etc. ). 
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After starting the data recording equipment, the next step of the completely automatic 
run sequence took electrical calibrations on all pressure transducers. Immediately fol- 
lowing this, the test tank was pressurized over a predetermined time period to the nom- 
inal operating pressure of 34. 47x10^ newtons per square meter (50 psia). Tank pressure 
was held constant for about 30 seconds to stabilize internal temperatures. The tank 
expulsion period was then started. Approximately 90 percent of the total tank volume was 
expelled at a constant volumetric flow rate. The expulsion period was stopped when a 
hot-wire liquid level sensor located at the 95-percent ullage level indicated gas. The last 
step of the automatic run sequence was the stopping of all data recording equipment. The 
test tank was then vented and refilled with liquid hydrogen for the next expulsion. 

Additional ramp pressurization runs, with no expulsion, were made for three dif- 
ferent tank ullage levels. The only deviation in the operating procedure for these runs 
was that the liquid outflow value was locked shut. As a result, after the ramp period, the 
hold period was simply extended as long as was desired. 


DATA REDUCTION 
Physical Description of Problem 


An initially vented tank containing a two -phase one -component cryogenic fluid was 
pressurized from 1 atmosphere to a new pressure by adding gas. The system was then 
allowed to reach a new equilibrium condition at which time liquid outflow was started. 
During this expulsion period, pressurant gas (at constant temperature) was added to the 
tank at a rate that maintained a constant tank pressure while expelling the liquid at a 
desired rate. The amount of pressurant gas used during the expulsion phase is dependent 
on (1) the volume of liquid displaced with no heat or mass transfer, (2) the heat transfer 
to the tank wall and liquid, and (3) the amount of mass condensed or evaporated. 

The main parameter used in the comparisons was the ratio of the ideal pressurant 
requirement to the actual pressurant requirement. The ideal pressurant wad determined 
under the assumption that the incoming pressurant gas did not exchange energy or mass 
with the surroundings. Under this assumption, the ideal pressurant required for the 
initial pressurization of the tank was determined by the relation 



MP qVq 

zrt g 



( 1 ) 


The ideal pressurant required for the expulsion period was determined by the relation 
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( 2 ) 


MPiVu 

Mj 

1 ZRT q 

Mass Balance 

A mass balance was performed on the ullage volume from an initial time tj to a 
final time tj as follows: 


M U,f “ M U,i + M G, i ■ 


O i Ml . £ 


(3) 


A discussion of how the terms of equation (3) were determined is now presented. 
Pressurant gas added (Mq j_f)- - The weight of the actual pressurant gas added 

from any initial time tj to any final time t f was determined by numerical integration of 
the gas orifice equation 


M, 


G,i-f 


f 


yD 2 C Vp Ap* dt 


(4) 


3 

Note that (tj - h ) is the time necessary to expel AVy (m ) of liquid. 

Ullage mass . - The initial ullage mass i and the final ullage mass M u f were 
obtained by numerical integration of the particular density profiles as follows: 


N4 


M tt . = / p dV p„ AV where p =f(T) 

U, i -Ar / ^ n n 

V U, i n=l 

r N f 

M U f = L P dV “ S p n AV n where p = f (T) 
’ V U,f n=l 


(5) 

( 6 ) 


The internal tank volume was considered as 57 (corresponding to thermopile location) 
horizontal disk segments. Each of these segments was in turn divided radially into a 
series of concentric rings, the number of which depended on the location of radial temper- 
ature sensors and the vertical position of the disk segment being considered. These 
rings (339 in all) comprised the V n ’s in the previous calculations. In this manner, 
vertical temperatures as well as radial temperature gradients could be incorporated into 
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the mass calculations. The position of the liquid level prior to and after expulsion deter- 
mined the number of gas volume rings (N^ and Nj) used in the ullage mass calculations. 
Mass transfer. - The mass transfer was calculated from equation (3) as a result of 


knowing My j, My and M^, that is, 


M t, i — f M U,i + M G, i— f " M U,f 

If M t . was a positive quantity, mass was considered leaving the ullage volume (e.g. , 
condensation). 


Energy Balance 

For the thermodynamic system consisting of the entire tank and its contents (tank + 
ullage gas + liquid), the first law of thermodynamics for an increment of time dt maybe 
written as 


y2 

dU T = (6 M g )(u g + P G V 0 + JL * Z G j - (6 M l )(u l + P L V L + _±i + Z L ) + «Q - 6W 


( 7 ) 


The kinetic and potential energy terms are small in comparison with the other energy 
terms and will be neglected in this development. If h = u + pv is substituted, equation (7) 
becomes 


dU T = (5M G )h Q - (6M L )h L + 6Q - 6W 


(8) 


For this system, there is no external work done on the system so 6W = 0. The final 
form of equation (7) is 

dU T = (6M Q )h G - (6M L )h L + 6Q (9) 

Equation (9) can be integrated over any time period. The physical interpretation of 
the quantities in equation (9) is as follows: 
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( 10 ) 



Change in 
system 
energy 
(Tank + gas 
+ liquid) 



Energy input 
by pressurant 
gas inflow 


Energy leaving 
through liquid 
outflow 


Energy from 
environment 
(Heat leak from 
conduction, 
convection, and 
radiation) 


The terms of equation (10) were evaluated in the following manner. 

Energy input b y pressurant gas inflow . - The first term in equation (10) may be 
evaluated as follows: 



M G h G 


dt ~ 


n=(t f -t./At) 



M G h G At 


(ID 


The pressurant flow rate was determined from equation (4). The specific 

enthalpy of the inlet gas was evaluated at the inlet temperature and pressure at each time 
increment At. 

Energy l eaving by liquid outflow . - 


/ 


M L h L dt =* 


n=(t f -t./At) 

% 


M L h L At 


( 12 ) 


The liquid flow rate was determined from the turbine flowmeter. The specific 
enthalpy of the liquid was evaluated at the outlet temperature. 

Energy input f rom environment. - The rate of energy input into the tank from the 
environment was assumed to be the same for all cases and was determined from a boiloff 

O 

test. This test indicated a nominal value of 0. 732x10 joules per second (0. 694 Btu/sec) 
should be used. This value which includes heat input by radiation, convection, and con- 
duction through pipes and supports was in all test cases less than 7. 0 percent of the 
energy added to the tank by the pressurant gas: 



dt 0. 732xl0 3 (t f 




(13) 
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Change in system ene rgy, - The change in system energy can be separated into the 
categories of (1) change in ullage energy, (2) change in liquid energy, and (3) change in 
the wall energy: 


dU T = dUjj + dU L + dU w 


( 14 ) 


Change in ullage energy . - The change in the ullage energy over any given time inter- 
val (tj — tj) is obtained by subtracting the internal energy of the ullage at time t. from 
the internal energy at time t^: 

/ U tf 

du u = ( u u) t - ( u u), <u> 

tf ti 

l i 

Making use of the relation U = H - PV gives 


f p( h -®) AV -Z 

Ju t v f \ P/ v t 



(16) 


The ullage gas density and enthalpy are functions of pressure and temperature. 
Therefore, by knowing the pressure and temperature profiles at times tj and t. one can 
evaluate the change in ullage energy. 

Change in liquid energy . - The change in energy of the liquid in the tank can be deter- 
mined in a manner similar to the change in ullage energy: 



Here again the liquid density and enthalpy are functions of pressure and temperature. 
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C hange in wall energy . - The change in wall energy was determined by applying the 
first law of thermodynamics to an element of the wall: 


/ 


' £dU w = AU w = AM w 


/ 


C v dT where C y = C y (T) 


The total change of the wall is then 

M„_ M 


w w /*T„ 

AU w, T = E ( AU w) n ' E AM W J C v (T) 


dT 


For convenience, equation (14) is substituted into equation (10): 

jf 1 ^ ( U u + U w + U l) dt - jf 1 “G h G dt - jf £ “L h L dt + J ' ‘ « 


Rearranging terms gives 


f 1 («G h G + «) dt = f 1 (“L h L dl + dU L ) + f 1 dU u + f 

«/t. v v j */t i t v J Jt J t. 


dU 


w 


Total energy 

Total change in 

Total 

Total change 

added 

liquid in tank 

change 

in wall 

(au t ) 

+ liquid 

in ullage 

energy 


expelled energy 

energy 

(AU W ) 


(aU L ) 

(AU n ) 


and dividing through by AUrp gives 


AUx AU tt AU w 

i = — — + n + w 

AU t AU t AU t 


(19) 


(20) 


( 21 ) 


( 22 ) 


(23) 


The data presented herein are in the form of these ratios. These ratios show the relative 
distribution of the total energy input. 
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TABLE I. - EXPERIMENTAL AND ANALYTICAL MASS BALANCE RESULTS 


Run 

Injector 

Tank 

Volume Inlet Ramp 

Hold Expulsion 

Tank 


Ramp period 



Expulsion period 



pressure, 

die- gas time, 

time, time, 

cycle 












N/m 2 

charge, temper- sec 

sec sec 

time, 

Initial 

Mass MasB 


Hold period 


Mass 

Mass Final 

Mass 

Mass 




3 _ i 



ullage 

added transfer 




added 

transfer ullage 

added 

transfer 




K 


sec 

mass, 

during during 

Ullage 

Mass Mass 

Ullage 

during 

during mass, 

during 

during 







kg 

ramp, ramp, 

mass 

added transfer 

mass 

expulsion, 

expulsion, kg 

expulsion, 

expulsion, 








kg kg 

after 

during during 

after 

kg 

kg 

kg 

kg 









ramp, 

hold, hold, 

, hold. 













kg 

kg kg 

kg 



(a) 

(a) 


14 

Hemisphere 

34. 7X10* 

1.50 

271 

15.0 

20.9 

261.3 

297. 2 

0. 167 

0.049 

0.029 

0. 188 

0.042 

0. 025 

0. 205±0. 002 

1. 369+0. 001 

-0. 248+0. 119 1. 822+0. 058 

1.308 

-0.178 

15 

[ 

34.7 

1.45 

266 

15.4 

20.3 

295.0 

330.7 

.224 

.054 

-.009 

.287 

.048 | 

.021 

. 314±0. 005 

1.382+0. 001 

-.119+0.074 1.816+0.062 

1.285 

-.202 

17 


34.9 

1.57 

271 

13.4 

22.3 

474.7 

510.4 

.079 

.031 

-.002 

.112 

.029 : 

.006 

. 134±0, 001 

1. 548+0. 003 

-.286+0.070 1.968+0.077 

1.545 

-.326 

18 


34.9 

1.56 

278 

13.2 

22.6 

204.1 

239.9 

.103 

.028 

.008 

.122 

.021 

-.001 

. 145+0. 001 

1.271+0.002 

-.153+0.045 1.569+0.041 

1. 234 

-.120 

20 


34.8 

1.58 

280 

12.7 

23.1 

132.4 

168.2 

.053 

.025 

.009 

.087 

.014 

-.010 

. 111±0. 001 

1. 121+0. 001 

-.109+0.036 1.342+0.026 

1.101 

-.083 

60 


35 7 

0 


32 0 

158 5 


190 5 

063 

031 

004 

QQO 

022 

- 228 

341±0 004 





63 


35.2 

1.59 



14.7 

22.4 

531.9 

569.4 

.061 

.034 

.027 

.068 

.027 

.003 

. 092+0. 001 

1. 565+0. 002 

-.355+0.09 1 2 . 012+0.082 

1.538 

-.382 

64 


34.7 

1.59 



14.2 

22.8 

377.6 

414.7 

.064 

.022 

.014 

.071 

.016 

-.025 

. 112+0. 001 

1. 410+0. 001 

-.255+0.067 1.777+0.063 

1.359 

-.248 

65 



1.58 



14.3 

22.6 

354.5 

391.4 

.079 

.019 

.006 

.093 

.014 

1 .009 

. 098±0. 001 

1.354+0. 001 

-.266+0.080 1.718+0.057 

1.312 

-.225 

67 



1.56 



15.6 

21.6 

263.7 

300.9 

.089 

.034 

.023 

.100 

.029 

' .010 

. 120+0. 001 

1.285+0. 001 

-.218+0.076 1.623+0.049 

1.212 

-.191 

68 



1.58 

1 


15.0 

22. 1 

149.0 

186.1 

.068 

.026 

.002 

.092 

.021 

.012 

. 102+0. 001 

1.057+0.002 

-.136+0.060 1.295+0.028 

1.011 

-.124 

69 


3F» 7 

0 

1 fi7 

32 R 

177 9 


210 7 

089 

039 

036 

093 

028 

- 041 

161+0 002 





70 


34.7 

0 



29 4 

150 1 


179 5 

073 

028 

- 015 

116 

019 

- 043 

176+0 002 





71 



1.58 



14.3 

23.6 

406.3 

442.2 

.071 

.039 

.036 

.077 

.027 

-.020 

. 123±0. 001 

1. 705+0. 003 

-.235+0.086 2.063+0.081 

1.679 

-.324 

73 



1.57 



13.8 

24.1 

288.9 

326.8 

. 100 

.035 

.001 

.135 

.023 

.014 

. 144±0, 001 

1. 652+0. 002 

-.129+0.068 1.926+0.060 

1.601 

-.166 

74 



1.58 



13.4 

24.5 

142.8 

180.7 

,061 

.033 

.0005 

.093 

.020 

-.004 

. 116±0. 001 

1. 426+0. 001 

.056+0.010 1.486+0.023 

1.452 

.061 

75 



0 

1 


31 0 

190 9 


221 9 

066 

D39 

007 

097 

023 

- 046 

169+0 002 


























484 

Radial 

34. 1X10 4 

1.54 

278 

14.3 

25.4 

336.5 

376.2 

0. 127 

0. 023 

0.006 

0. 143 

0. 022 

0 

0. 165±0. 006 

1. 402+0. 001 

-0. 274+0. 037 1. 841+0. 070 

1.372 

-0.256 

486 



1.55 

288 

13.5 

26.3 

277.3 

317.1 

.127 

.018 

.012 

.134 

.013 

-.026 

. 173+0. 002 

1.344+0.001 

-.196+0.027 1.714+0.061 

1.320 

-.207 

487 



1.54 

291 

14.3 

25.6 

191. 1 

231.0 

.133 

.024 

.005 

.161 

.019 

.015 

. 166+0. 006 

1.221+0.001 

-. 150+0. 033 1. 536+0. 040 

1.191 

-.146 

488 



1.54 

293 

14.5 

25.1 

134.8 

173.9 

.127 

.024 

.008 

.143 

.020 

.001 

. 162 ±0. 005 

1. 088+0. 001 

-. 118+0. 026 1. 369+0. 027 

1.091 

-.107 

489 



1.54 

284 

14.0 

25.4 

402.9 

442.3 

.120 

.022 

-.003 

.145 

.019 

o 

o 

. 167±0. 007 

1.435+0.002 

-.302+0.067 1.904+0.078 

1.422 

-.303 

490 

Straight 

34. 3xl0 4 

1.55 

290 

15.4 

24.4 

402.5 

442.3 

0. 121 

0. 028 

0.021 

0. 128 

0. 021 

-0. 027 

0. 176 ±0, 002 

1. 292+0. 002 

-0.329+0.033 1.798+0.059 

1.279 

-0.283 

491 

pipe 


1.55 

289 

14.7 

24.7 

334.7 

374.1 

. 114 

.026 

.007 

.133 

.019 

. 005 

. 147±0. 002 

1.226+0. 002 

-.349+0.047 1.722+0.068 

1.231 

-.231 

492 



1.54 

289 

14.7 

24.7 

275.5 

314.9 

.139 

.024 

.027 

.137 

.020 

’ -.020 

. 176+0. 002 

1.156+0. 001 

-.234+0.026 1.566+0.045 

1. 149 

-.185 

493 



1.54 

289 

14.7 

24.8 

190.3 

229.8 

.142 

.027 

.025 

.143 

.020 

| -.001 

. 164+0. 002 

1.014+0. 001 

-.223+0.019 1.400+0.031 

1.054 

-. 126 

494 



1.55 

291 

14.8 

24.6 

135.6 

175.0 

.119 

.026 

.013 

.133 

.019 

.005 

. 147+0. 002 

. 895+0. 001 

187+0. 014 1. 229+0. 022 

.964 

-.088 

1 

Hemisphere 

34.2x10* 

1.51 

294 

23,2 

29.3 

394.0 

446.5 

0. 112 

0. 033 

0.027 

0.114 

0. 007 

-0. 039 

0. 165+0. 003 

1. 338+0. 001 

-0.196+0.049 1.700+0.064 

1.256 

-0.275 

2 


34.2 

1.52 

294 

22.7 

29.5 

332.5 

384.7 

.092 

.033 

.014 

.111 

.013 

-.015 

. 139+0. 001 

1.310+0. 001 

-.162+0.031 1.611+0.049 

1.211 

-.253 

3 


34.1 

1.52 

285 

23,0 

29.6 

275.9 

328.5 

.094 

.039 

.019 

.114 

.015 

0 

. 129+0. 001 

1. 272+0. 001 

-.150+0.032 1.550+0.047 

1.175 

-.237 

4 


34.2 

1.54 

287 

22.9 

29.5 

193.9 

246.3 

.078 

.034 

-.006 

.117 

.014 

.014 

. 140+0. 001 

1. 192+0. 002 

-.110+0.032 1.420+0.036 

1.116 

-.156 

5 


34.1 

1.53 

2 

189 

22.8 

29.7 

138.2 

190.7 

.082 

.032 

.015 

.099 

.012 

, .033 

. 144+0. 001 

1.065+0. 001 

. -.068+0.018 1.279+0.025 

1.066 

-.093 


l Values predicted by analytical program. 



TABLE H. - EXPERIMENTAL AND ANALYTICAL ENERGY BALANCE RESULTS 


Run 


Injector 


Inlet 


Expulsion 


Expulsion period 



gas 

temper- 

ature, 

K 

time, 

sec 

Energy added 
by pressurant 
gas, 

J 

Energy 
[ added by 
environment, 

j 

Energy gained 
by tank wall, 

AU w- 

J 

Energy 

gained by 

ullage, 

AU 

u 

Energy gained 
by liquid, 

au l> 

J 

Experimental 

Analytical 

Experimental 

Analytical 

Experimental 

Analytical 

14 

Hemisphere 

271 

261.3 

530.2±0.4xl0 4 

506. 3xl0 4 

19. lxlO 4 

268. 6xl0 4 

263. 7xl0 4 

169. 1±2. lxlO 4 

109. 0±19. 7xl0 4 

72. 8xl0 4 

15 



266 

295.0 

529. 8±0. 5 

492.2 

21.6 

274.6 

266.9 

159. 9±1. 9 

91. 0±22. 6 

85.9 

17 



271 

474.7 

599. 5±1. 1 

598.3 

34.8 

313.5 

326.8 

182. 8±2. 3 

107. 7 ±26. 6 

138.6 

18 



278 

204.1 

501. 8±0. 6 

486.7 

14.9 

234.9 

249.0 

163. 6±1. 9 

78. 7±17. 2 

57.8 

20 



, 280 

132.4 

445. 9±0. 4 

437.9 

9.7 

189.0 

202.8 

156. 6±1. 7 

73. 1±13. 7 

39.0 

63 



333 

531.9 

688. 7±0. 8 

676.3 

38.9 

413.9 

481.4 

189. 7±2. 4 

123. 8±15. 3 

162.8 

64 





377.6 

685. 5±0. 6 

660.1 

27.7 

409.5 

443.4 

178. 0±2. 3 

125. 0±13.2 

114.8 

65 





354.5 

689. 0±0. 6 

667.0 

26.0 

402.3 

443.4 

175. 2±2. 1 

123.0±12.6 

107.6 

67 





263.7 

623. 6±0. 4 

588.1 

19.3 

378.1 

372.4 

168. 5 ±2. 1 

122. 1±10. 7 

80.4 

68 





149.0 

544. 0±0. 9 

520.1 

10.9 

251.3 

290.8 

156. 5±1. 8 

109. 6±15. 5 

45.4 

71 



167 

406.3 

449. 4±0. 7 

442.2 

29.8 

194.6 

181.2 

186. 9±2. 4 

111. 4±12. 7 

123.3 

73 



167 

288.9 

413. 7±0, 4 

400.9 

21.2 

162.2 

143.0 

178. 2 ±2. 2 

80. 5±11. 2 

88.3 

74 

i 




167 

142.8 

347. 0±0. 3 

353.2 

10.5 

108.2 

101.4 

159. 8±1. 7 

90. 7±14. 9 

43.3 

484 

Radial 

278 

336.5 

550. 6±0. 5xl0 4 

538. 5xl0 4 

24. 6xl0 4 

283.7X10 4 

295. OxlO 4 

172. 1±2. 2xl0 4 

97. 5±22.2xl0 4 

102. 7x1 0 4 

486 



288 

277.3 

539. 5±0. 4 

529.8 

20.3 

263.9 

286.6 

166. 0±2. 0 

73. 7 ±20. 7 

84.9 

487 



291 

191.1 

496. 9±Q. 4 

484.5 

13.4 

222.2 

249.3 

157. 8±1. 8 

63. 9±17. 0 

65.0 

488 



293 

134.8 

445. 2±0. 3 

446.1 

9.9 

154.7 

208.4 

150. 3±1. 7 

55. 0±14. 9 

41.4 

489 



284 

402.9 

569, 3±0. 6 

563.6 

29.5 

309.9 

316.6 

174. 9±2. 3 

104. 7 ±2 4. 4 

122.5 

490 

Straight 

290 

402.5 

521. 6±0. 4xl0 4 

515.9X10 4 

29. 5x1 0 4 

257.0X10 4 

297. lxlO 4 

169. 9±2. lxlO 4 

77. 8±24. 7xl0 4 

122. OxlO 4 

491 

pipe 

289 

334.7 

494. 2±0. 6 

496. 2 

24.5 

233.6 

271.3 

167. 2 ±2. 2 

69. 6±22. 4 

101.8 

492 



289 

275.5 

485. 1±0. 6 

481.7 

20.2 

210.5 

251.5 

158. 7±1. 7 

73. 3±20. 2 

83.7 

493 



289 

190.3 

410. 5±0. 5 

426.5 

13.9 

168.1 

208.7 

151. 3±1. 5 

54. 4±17. 2 

57.9 

494 



291 

135.6 

365. 3±0. 3 

393.4 

9.9 

130.4 

178.4 

145. 2±1. 6 

48. 3±14. 8 

41.3 

1 

Hemisphere 

294 

394.0 

546. 7±1. 5x10 4 | 

512.8X10 4 

28. 8xl0 4 

294. 8xl0 4 

269. lxlO 4 

169. 9±2.2xl0 4 

76. 1±24. 4xl0 4 

126. 5xl0 4 

2 



294 

332.5 

529. 5±1. 4 

489.3 

24.3 

273.6 

247.5 

167. 7±2. 0 

80. 4±22. 2 

105.5 

3 



285 

275.9 

514. 0±2. 1 

474.4 

20.2 

251.8 

231.5 

165. 7±2. 0 

85. 4±20. 2 

87.9 

4 



287 

193.9 

484. 1±0. 8 

453.1 

14.2 

225.8 

205.2 

160. 6±1. 9 

83. 0±17. 2 

62.3 

5 



289 

138.2 

434. 4±5. 5 

434.4 

10.1 

197.6 

171.4 

152. 1±1. 7 

65. 6±14. 9 

44.5 


Error Analysis 


An analysis was performed to determine the magnitude of probable error which could 
be present in the integrations of equations (4) to (6) and in the determination of all terms 
in the energy balance equation (13). Probable error is defined as follows: There is a 
50-percent probability that the error will be no larger than the value stated. This analy- 
sis considered the errors introduced by the inaccuracies of temperature transducers as 
well as the tank pressure sensor. These calculations were performed for all runs for the 
expulsion period. The results of this analysis are included with the tabular data in 
tables I and n. 


RESULTS AND DISCUSSION 

The main criterion used to compare the effectiveness of the various operating param- 
eters on the amount of pressurant gas used is the nondimen sional ratio Mj/Mq, where 
Mj is the ideal pressurant mass required to expel a given volume of liquid at a given 
inlet gas temperature and tank pressure. The actual pressurant requirement for the 
same conditions is Mq. A higher Mj/Mq ratio means less energy and mass exchange. 

It does not necessarily mean a lower absolute pressurant requirement Mq as will be 
illustrated later in this section. 

For no environmental heating, a value of Mj/Mq equal to unity implies that no heat 
or mass transfer occurs. Therefore, from equation (23), the AUjj/AILj, ratio is also 
equal to unity; that is, all the energy (AU T ) added to the tank during expulsion appears as 
an increase in ullage energy (AUjj). Accordingly, any value of Mj/Mq or AUjj/aU t 
less than unity means energy is lost by the ullage system. This loss of ullage energy 
would then appear as an increase in tank wall heating and/ or liquid heating AU^/ AUrp 
and AU l /AU t . 

The discussion of results will first present the effects of the various operating param- 
eters on the Mj/Mq ratio, for the expulsion period only, which was of primary interest, 
followed by the mass transfer M^/Mq results. Then the results of the energy balances 
will be presented in an attempt to point out the major reason for the Mj/Mq or AUjj/AU T 
ratio being less than one. Finally a comparison will be made between the experimental 
results and the analytically predicted results to determine the validity of the analytical 
program. The analytical results are presented in the figures along with the correspond- 
ing experimental results. The comparison between experimental and analytical results 
will be given in terms of an average deviation. Average deviation is defined as 
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I 



TABLE m. - DEVIATIONS BETWEEN EXPERIMENTAL ANALYTICAL RESULTS 


Run 

Injector 

Inlet gas 

Tank wall 

Percent deviation between experimental 




temperature, 

thickness, 


and analytical results 3. 





K 

cm 

Mj/M g 

M i/ M GN 

au w /au t 

au l /au t 

14 

Hemisphere 

271 

0.762 

+8.09 

-8.50 

+1. 84 

+28.6 

15 



266 



+7.18 

-.90 

+3.01 

+1.27 

17 



271 



+10. 45 

-6.20 

-4.25 

-28.80 

18 



278 



+6.33 

-4.90 

-6.37 

+20. 40 

20 



280 



+5.86 

-4. 00 

-8.00 

+ 46.60 








Av +7.58 

Av -4.90 

Av 4. 70 

Av 25. 13 

63 

Hemisphere 

-333 

0.762 

+18.69 

0 

-16. 16 

-30.60 

64 







+12.45 

-3.50 

-8.04 

+2.86 

65 







+11.97 

-5.50 

-8.06 

+9.88 

67 







+8.36 

-7.00 

+1.36 

<*33. 70 

68 







+6.86 

-5.30 

-15.66 

+57.10 








Av +11.67 

Av -4.30 

Av 9.86 

Av 26.83 

71 

Hemisphere 

167 

0.762 

+21.28 

+10.60 

+ 6.90 

-9.92 

73 







+18.02 

+11.60 

+11.80 

-10. 80 

74 







+16.69 

+21.00 

+ 6.27 

+52.40 








Av +18. 66 

Av +14. 40 

Av 8.32 

Av 24. 4 

484 

Radial 

278 

0. 762 

+13.67 

-2.90 

-3.85 

-5.31 

486 



288 



+11.78 

-1.00 

-8.03 

-16.03 

487 



291 



+8.86 

-2. 50 

-12.41 

-1.60 

488 



293 



+9.20 

-.60 

-11.19 

-+1. 80 

489 



284 



+16.68 

-.40 

-2.32 

-17.14 








Av +12. 04 

Av -1.50 

Av 7.56 

Av 12.98 

490 

Straight 

290 

0.762 

+22.20 

+3.62 

-28. 53 

-58.90 

491 

pipe 

289 



+21.30 

-1.05 

-19.69 

-23.90 

492 



289 



+20.30 

+4.37 

-19.42 

-10. 00 

493 



289 



+20.50 

+3.62 

-10. 88 

-8.60 

494 



291 

i 


+21. 50 

+5. 16 

-7.93 

+13.95 








Av +21. 16 

Av +3. 56 

Av. 17.29 

Av 23. 07 

1 

Hemisphere 

294 

0 . 

409 

+11.56 

0 

+9.00 

-66.70 

2 



294 



+10. 19 

-1.00 

+9.52 

-31.70 

3 



285 



+9.89 

-.90 

+8.07 

-2.50 

4 



287 



+6.44 

-2.30 

+9.28 

+25.15 

5 


r 

289 



+8.06 

-2.20 

+13.06 

+32.40 








Av +9. 23 

Av -1.30 

Av 9.78 

Av 31.69 


a Analysis underpredicted, +; analysis overpredicted, 
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[El 

/Experimental^ /Analytical \ 
y ratio / \ ratio / 

1 

^Experimental^ 
y ratio ) 

1 

N 


( 100 ) 


(24) 


where N is the number of data points in a given set of operating conditions (i. e. , when 
using a hemisphere injector for a constant inlet gas temperature of 167 K (301° R), N 
would be 3 for the data presented in fig. 8). For convenience, all deviations between the 
experimental and analytical results are summarized in table in. 

The operating parameters (e.g., inlet gas temperature, injector geometry, tank 
wall thickness) and the major experimental and analytical results are summarized in 
tables I and II. Table I gives experimental and analytical mass balance results, while 
table II gives the corresponding energy balance results. 

Temperature effect on pressurant gas requirements using a diffusing injector (hemi- 
sphere). - The effect of inlet gas temperature is shown in figure 8 on the basis of Mj/Mq 
for various expulsion times. Expulsion time is the total time required to expel liquid 
from a 5 -percent ullage to a 9 5 -percent ullage. Therefore, each data point represents a 
complete expulsion. As may be seen in the figure, the Mj/M g ratio decreases for 
increasing inlet gas temperature. This implies a larger percentage of energy that is con- 
tained in the pressurant gas is lost to the tank wall and liquid as the inlet gas temperature 
is increased. The average decrease in Mj/Mq for all expulsion times in going from a 
167 K (301° R) inlet gas temperature to a 273 K (491° R) inlet temperature is 27. 0 percent. 



Figure 8. - Comparison of ideal pressurant requirement to actual pressurant requirement ratio as 
function of expulsion time for three inlet gas temperatures. Hemisphere injector; wall thickness, 
0. 762 centimeter (0. 30 in. ); tank pressure, 34. 47X10 4 newtons per square meter (50 psia). Curves 
faired through experimental data. 
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Figure 9. - Comparison of actual pressurant requirements for various expulsion times for three inlet 
gas temperatures. Hemisphere injector; tank wall thickness, 0.762 centimeter (0.30 in. ); tank 
pressure, 34.47xl(r newtons per square meter (50 psia). Curves faired through experimental data. 


A comparison of the actual pressurant requirements Mq for the three inlet gas temper- 
atures for various expulsion times is shown in figure 9. For convenience, the various 
expulsion times (or constant liquid outflow rates) are also shown in terms of an equivalent 
thrust level based on the specific impulse of the RL10A3-3 engine (Ig = 444 sec). The 
actual pressurant requirement Mq decrease for increasing inlet gas temperature and 
increase for increasing expulsion times. The data show an average reduction of 1. 40 per- 
cent in pressurant gas requirements for a 10 K (18° R) increase in inlet gas temperature 
over the range of expulsion times. Although the absolute value of Mq decreases for 
increasing inlet gas temperatures, the Mj/Mq ratio also decreases because of the even 
greater decrease in ideal requirements. The average decrease in Mj/Mq in going from 
a 273 to 334 K (491° to 601° R) inlet gas temperature is 17. 20 percent. For a given inlet 
gas temperature, there is an increasing pressurant requirement (decreasing Mj/Mq) for 
increasing expulsion times. The longer the pressurant (ullage) gas is exposed to cold 
surroundings, the greater the loss in pressurant energy. 

The shaded symbols shown in figures 8 and 9 are the results as predicted by the 


25 



analytical program. The analysis underpredicts the ratio Mj/Mq, overpredicts the 
actual pressurant requirements, by an average of 12. 64 percent for the three inlet gas 
temperatures (see table m, column (a), for actual deviations). The scatter of the ana- 
lytical results is due to variations in the actual initial conditions (i. e. , initial ullage vol- 
ume, initial gas and wall temperatures, inlet gas temperatures, and pressure) that are 
used as input to the analytical program. 

It should be kept in mind that this modified analysis still neglects mass transfer. 
Therefore all the predicted pressurant requirements must be supplied from a source 
external to the system (through the inlet injector), whereas in the true physical case (as 
will be shown later) some of the mass added to the ullage is supplied by the evaporation of 
liquid. At best the analytically predicted pressurant requirement should be greater than 
the experimental value by an amount equal to the mass evaporated (or less than the exper- 
imental value by an amount equal to the mass condensed). 

If the comparison just presented were based on the net experimental mass added to 
the ullage (i. e. , actual pressurant through injector plus any mass evaporated), the 

analysis overpredicts the ratio by an average of 4. 60 percent for the 273 and 

334 K (491° and 601° R) inlet temperatures and underpredicts the Mj/Mq^- ratio by an 
average of 14. 40 percent for the 167 K (301° R) inlet gas temperature (see table in, 
column (b)). 

Because of this limitation imposed on the pressurant requirements as predicted by 
the analysis, the discussion of the comparisons between experimental and analytical pres- 
surant requirements will include results based on both the actual pressurant added through 
the injector Mq and the net mass added to the ullage Mq N . The later comparison will 
do more justice to the analytical program and its assumptions. 

A comparison of the ratio M^/Mq of mass transferred during expulsion to the actual 
pressurant added to the tank is presented in figure 10 for different expulsion times and 
the three inlet gas temperatures. In all but one run, the net mass transfer was evapor- 
ation and represented between 8 and 23 percent of the pressurant mass added to the tank 
during expulsion. For two of the three inlet gas temperatures, the M^/M^ ratio 
increases for increasing expulsion time. The same trend is present for the 167 K (301° R) 
temperature for expulsion times greater than 180 seconds. For a given expulsion time, 
the ratio increases for increasing inlet gas temperature. 

The AU n /AU T ratios of the energy increase in the ullage over the expulsion period 
to the total energy added to the tank is compared in figure 11 for different expulsion times 
and for the three inlet gas temperatures. For all runs, between 25 and 45 percent of the 
total energy that was added to the tank remains in the ullage after expulsion. For any 
given expulsion time, the AU' U /AU' T ratio decreases (indicating increased energy loss to 
the wall and liquid) for increasing inlet gas temperatures. Also, for a given inlet gas 
temperature, the AU u /AU t ratio decreases with increasing expulsion time. It should 
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Runs 



Figure 10. - Comparison of mass transfer to pressurant added ratio as function of expulsion 
time for three inlet gas temperatures. Hemisphere injector; tank wall thickness, 0. 762 centi- 
meter (0. 30 in. ); tank pressure, 34.47X10 4 newtons per square meter (50 psia). Curves faired 
through experimental data. 


Inlet gas Runs 

temperature, 

K(°R) 

□ 167 (301) 71, 73, 74 

O 273 (491) 14, 15, 17, 18, 20 



Total expulsion time, sec 

Figure 11. - Comparison of energy increase in ullage to total energy added ratio as function of expul- 
sion time for three inlet gas temperatures. Hemisphere injector; tank wall thickness, 0. 762 
centimeter (0. 30 in. ); tank pressure, 34.47X10 4 newtons per square meter (50 psia). 


be noted (table II) that the absolute value of AUy does not change significantly for the 
various operating parameters used during testing. The two terms used in the determin- 
ation of AU|j are ullage gas mass and specific internal energy. As the nominal value of 
inlet gas temperature is increased, the specific internal energy of the ullage gas neces- 
sarily is higher; however, the mass of gas in the ullage is significantly reduced. The net 
result of these two countering trends is to produce only a slightly decreasing value of 
AUjj for increasing inlet gas temperature. The mean increase in ullage energy for all 
runs was 166. 2xl0 4 joules (1576 Btu) with a standard deviation of llxlO 4 joules (104. 5 
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Figure 12. - Comparison of energy gained by wall to total energy added ratio as function of expulsion 
time for three inlet gas temperatures. Hemisphere injector; tank wall thickness, 0. 762 centimeter 
(0. 30 in. ); tank pressure, 34. 47x10^ newtons per square meter (50 psia). Curves faired through 
experimental data. 


Btu). Any trends in the AU^/ AU T ratio, therefore, depend mainly on variations in the 
total energy added (AU T ) to the tank due to variations in energy losses to the tank wall 
and liquid. 

The ratio AU W / AUrp of energy gained by the tank wall to the total energy added is 
compared in figure 12 for three inlet gas temperatures for different expulsion times. In 
general, between 30 and 58. 8 percent of the total energy added to the tank was gained by 
the tank wall over the range of conditions. As can be seen in table II and in figure 12, 
both the absolute value of AU W and the AU^y/AUrp ratio increase with increasing inlet 
gas temperature. The increase in AU W is due to a larger driving potential (AT) for heat 
transfer between the ullage gas and the tank wall. The total energy added to the tank 
AU t does not increase in the same proportion as AU W , which results in the increasing 
AU W /AU T ratio. Both the absolute value AU W and the AU W /AU T ratio also 
increase with increasing expulsion time. The increase in AU W for increasing expulsion 
time is due to the longer exposure of the tank wall to a heat source. 

Figure 13 is a comparison of the ratio of energy gained by the liquid to the total 
energy added to the tank AUp/ AUrp for different expulsion times for the three inlet gas 
temperatures. In all cases, between 15 and 25 percent of the total energy added to the 
tank appears as an increase in energy of the liquid. There are no identifiable trends in 
the experimental data with either increasing inlet gas temperature or increasing expulsion 
time. 
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Figure 13. - Comparison of energy gained by liquid to total energy added ratio as function of expul- 
sion time for three inlet gas temperatures. Hemisphere injector; tank wall thickness, 0. 762 centi- 
meter (0.30 in.) ; tank pressure, 34.47x10^ newtons per square meter (50 psia). 


The analytical predictions of the ratio AU^/AUrp are also presented in figure 13. 

(See table in for average deviations. ) The large discrepancy between analysis and exper- 
iment could be the result of the error in experimentally determining AU^. The probable 
error associated with each experimental determination of AU^ is between 8. 8 and 25 
percent (see table II). The analysis generally underpredicts the ratio below 340 to 380 
second expulsion and overpredicts the ratio above 340 to 380 second expulsion. While the 
actual data showed no trends with either expulsion time or inlet gas temperature, the 
analysis does. The analysis predicts decreasing AUjVAU T with increasing inlet gas 
temperature for constant expulsion time and predicts increasing AUp/AUrp with increas- 
ing expulsion time for constant inlet temperature. 

The results from the preceding discussion point out that in all cases between 73. 5 
and 85 percent of the total energy that was added to the tank is either absorbed by the 
tank wall or remains in the ullage. The correlation between the analysis and experimental 
data, therefore, depends largely on the ability of the analysis to predict final wall and 
ullage gas temperature profiles. These temperature profiles are, in turn, used to deter- 
mine the increase in wall and ullage energy and the final ullage mass. The ability to 
predict these temperatures explains the good agreement between experimental data and 
analysis reported in reference 1. 

Figures 14 to 16 present comparisons of experimental and analytical wall and ullage 
gas temperature profiles. Three particular runs were selected to represent cases where 

(1) the agreement between analytical and experimental values of AU^/AUrp was the best, 

(2) the analysis overpredicts AU w /AU t , and (3) the analysis underpredicts AU w /AU t . 

Figure 14 presents a comparison of experimental and analytical wall and ullage gas 
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Figure 14. - Comparison of analytical and experimental gas and wall temperatures at end of 261-second expul- 
sion. Hemisphere injector; inlet gas temperature, 271 K (488° R); tank wall thickness, 0. 762 centimeter 
(0.30 in.); tank pressure, 34.47x10^ newtons per square meter (50 psia) ; run 14. 


temperature profiles for a run where the ratios Mj/Mq and AU^/AlLp that were pre- 
dicted by the analysis were within 8. 09 and 1. 84 percent of the experimental results. The 
experimental gas temperatures shown on figure 14 are obtained from the vertical rake and 
indicate average radial temperatures at their respective vertical positions. In the 
absence of any mass transfer, the pressurant mass required for an expulsion could be 
determined as the difference between the final mass in the ullage and the initial mass 
prior to expulsion. Since one of the initial conditions for the analytical program is the 
initial experimental temperature profile (and, therefore, the initial ullage mass) prior to 
expulsion, the deviation in pressurant requirements would largely be the result of the pre- 
dicted final temperature profile. As seen in figure 14, the calculated gas temperatures 
are slightly lower than the experimental temperatures in the lower 18 percent of the 
ullage (where 35 percent of the experimentally determined ullage mass is concentrated) 
but much higher in the remaining ullage. However, due to the temperature -density 
relation, the mass present in the final ullage volume for the two profiles differ by only 
10 percent. The 10 percent deviation in final ullage mass plus the fact that the analysis 
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Figure 15. - Comparison of analytical and experimental gas and wall temperatures at end of 532-second expul- 
sion. Hemisphere injector; inlet gas temperature of 334 K (601° R); tank wall thickness, 0.762 centimeter 
(0. 30 in.); tank pressure, 34.47X10 4 newtons per square meter (50 psia); run 63. 


neglects mass transfer account for the 8. 09 percent deviation in ratio (8. 5 per- 

cent in Mj/m gn ratio). 

The analytically predicted wall temperatures (fig. 14) are very close to those 
observed experimentally in the lower half of the tank ullage but begin to get higher toward 
the top of the tank. Even though the analytical wall temperatures are higher in the upper 
portion of the tank wall, the analysis underpredicts AU^ by 1. 84 percent. All values of 
energy gained by the tank wall (AU^) include the energy gained by the tank lid. The tank 
lid is approximately 22 percent of the total tank mass. The final lid temperature, both 
predicted and actual are shown as the uppermost wall temperature (point zero in fig. 14). 
The predicted lid temperature is 15. 5 K (28° R) lower than the observed temperature. 

The underprediction of temperature for a large mass concentration together with the fact 
that the specific heat of the wall is a strong function of temperature apparently offsets the 
increase in wall energy due to overpredicting the temperatures in the upper portion of the 
tank wall. 

Figure 15 is a comparison of experimental and analytical wall and ullage gas temper - 
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Figure 16. - Comparison of analytical and experimental gas and wall temperatures at end 
of 289 -second expulsion. Hemisphere injector; inlet gas temperature, 167 K (301° R); 
tank wall thickness, 0.762 centimeter (0.30 in.); tank pressure, 34. 47X10 4 newtons 
per meter (50 psia); run 73. 


ature profiles for a run which had an inlet gas temperature of 334 K (601° R). Here again, 
the calculated ullage gas temperatures are lower than the experimental temperatures in 
the lower 25 percent of the ullage (where 46 percent of the ullage mass is concentrated) 
but much higher in the remaining ullage. The total final ullage mass is approximately the 
same for both profiles. As a result, the deviation between the analytical and experimental 
values of Mj/Mq is entirely due to the analysis neglecting mass transfer. 

The analytically predicted wall temperatures agree fairly well with the experimental 
temperatures except in the upper 25 percent of the wall where the analysis overpredicts 
temperatures as much as 58 K (104° R). For this run, the analysis overpredicts the 
energy gained by the wall. In this case, the underprediction of lid temperature (point 
zero in fig. 15) does not offset the increase in wall energy due to the large overprediction 
of wall temperatures in the upper portion of the tank wall. 

A comparison of experimental and analytical wall and ullage gas temperatures for a 
run which had an inlet gas temperature of approximately 167 K (301° R) is shown in fig- 
ure 16. Here the analytical wall and gas temperatures are very close to those obtained 
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However, the analytical prediction of the mass ratios Mj/Mq 


or 


experimentally, 

Mj/Mgn was poor (typical of the 167 K (301° R) prediction). The analysis underpredicts 
the energy gained by the wall as a result of underpredicting the temperature rise of the 
lid. 


From the preceding discussion it is evident that the analysis predicted accurately 
the pressurant requirements even though the analysis, in most cases, did not accurately 
predict the temperature profile in the ullage after the expulsion. The fair agreement 
between the analysis and observed requirements is because the "mass weighted average 



is approximately equal for both the predicted 


and observed final temperature profiles. 

Injector e ffect on pressurant requirements. - The effect of injector geometry on the 
Mj/Mq ratio for different expulsion times for an inlet gas temperature of 284±6. 0 K 
(511±11° R) is shown in figure 17. The hemisphere and radial injectors (diffuser types) 
have approximately the same Mj/Mq ratio for the various expulsion times. The 
straight pipe, however, has lower pressurant requirements (a greater Mj/Mq ratio) 
than the diffuser -type injector for all expulsion times. The straight pipe injector required 
20 percent less pressurant for a 135-second expulsion and 8. 9 percent less for a 402- 
second expulsion. At some expulsion time greater than 500 seconds, the straight pipe and 
diffuser type injectors would appear to have the same pressurant requirement. At this 
point, the effect of the injection pattern of the straight pipe would be minimized due to the 
low injection velocities. It should be noted here that these are approximately the same 
results that were obtained in reference 3 with a small cylindrical tank. The reasons for 
the decreased pressurant requirements (increased Mj/Mq ratio) when using the straight 
pipe injector are (1) less heat transfer to the tank walls and (2) a greater evaporation of 
liquid. These two facts will be discussed further later in this section. All three injectors 
have increasing pressurant requirements (decreasing Mj/Mq ratio) for increasing expul- 
sion times. 

The analytical predictions for the three injectors can be approximated by a single 
curve as shown in figure 17. Since inlet gas temperature and tank pressure as functions 
of time and liquid flow rate are the only required inputs to the analytical program, it is 
not surprising that the Mj/Mq ratio is not a function of injector design. The analytical 
predictions of Mj/Mq are better for the diffuser type injectors because this type best 
approximates the one -dimensional flow assumption of the analysis. The average deviation 
between the analytical and experimental values of Mj/Mq for the diffuser type injectors 
is 9. 81 percent compared to 21. 16 percent for the straight pipe injector. 

The agreement between analysis and experimental values of m j/ m gn are good for 
all three injectors. The average deviation for all three injectors is 3. 67 percent. The 
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Figure 17. - Comparison of ideal pressurant requirement to actual pressurant require- 
ment ratio as function of expulsion time for three injector geometries. Inlet gas tem- 
perature, 284±6 K (511°±11° R) ; tank wall thickness, 0. 762 centimeter (0. 30 in. ); tank 
pressure, 34.47x10* newtons per square meter (50 psia). 


improved agreement, especially for the straight pipe injector, is probably the result of 
compensating limitations or assumptions of the analytical program. This compensation 
will be discussed later in this section. 

The comparison of the M^/M^ ratios for different expulsion times for the three 
injector geometries is shown in figure 18. In general, all injectors show increased 
M t/ M G ( increased evaporation) with increasing expulsion time. The straight pipe injec- 
tor has a higher M^/Mq ratio (more evaporation) for all expulsion times than the dif- 
fuser type (93. 5 percent higher than the radial injector at a 135-sec expulsion and 21. 4 



Total expulsion time, sec 

Figure 18. - Comparison of a mass transfer to pressurant added ratio as function of 
expulsion time for three injector geometries. Inlet gas temperature, 284±6 K 
(511°+11° R); tank wall thickness, 0. 762 centimeter (0. 30 in. ); tank pressure, 
34.47X10 4 newtons per square meter (50 psia). 
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Figure 19. - Comparison of energy increase in ullage to total energy added ratio as func- 
tion of expulsion time for three injector geometries. Inlet gas temperature 284+6. OK 
<511°±11° R); tank wall thickness, 0. 762 centimeter (0. 30 in. ); tank pressure, 34.47X10 4 
newtons per square meter (50 psia). 


percent higher at a 402 -sec expulsion). Of the diffuser type injectors, the radial injector 
has the higher M^/M^ (11. 3 percent higher than the hemisphere at 135- sec expulsion 
and 23. 5 percent higher at 402-sec expulsion). 

The comparison of the AU^/AU^, ratio for different expulsion times for the three 
injector geometries is shown in figure 19. In general, AU^/ AU^, decreases with 
increasing expulsion time for the three injector geometries. The diffuser type injectors 
(hemisphere and radial) have approximately the same AU u /aU t ratio for the various 
expulsion times. Figure 19 indicates that, for the diffuser type injectors, between 29 
and 34 percent of the total energy added to the tank remains in the ullage at the end of the 
expulsion periods. The straight pipe has an average 9. 3 percent higher AU u /AU t ratio 
than the diffuser type injectors. For the straight pipe, between 31. 0 and 39. 0 percent of 
the total energy added to the tank remains in the ullage after expulsion. 

The ratio of energy gained by the tank wall to the total energy added to the 0. 762 



Figure 20. - Comparison of energy gained by wall to total energy added ratio as function of 
expulsion time for three injector geometries. Inlet gas temperature, 284+6 K 
(5U°±H 0 R); tank wall thickness, 0.762 centimeter (0.30 in. ) ; tank pressure, 

34. 47x10^ newtons per square meter (50 psia). Curves faired through experimental data. 
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centimeter (0. 30 in. ) thick tank wall AU W /AU T for different expulsion times for the 
three injector geometries is shown in figure 20. The two diffuser type injectors have 
approximately the same AU w /AU t ratios for the various expulsion times. The data 
presented can be represented by a single curve. The straight pipe injector causes less 
wall heating and a lower AU w /AU t ratio than the diffuser types - 15. 60 percent lower 
at a 135-second expulsion and 9. 90 percent lower at a 402-second expulsion. All three 
injectors have increasing wall heating and increasing AU w /AUrp for increasing expulsion 
times. The increased wall heating is because the wall is exposed to a heat source for a 
longer duration for the longer expulsion times. The reduced wall heating when using the 
straight pipe is due to lower, nearly constant, axial temperature gradients in the upper 
portion of the ullage together with rather large radial temperature gradients that were 
lowest near the tank wall and highest at the center of the tank. 

In general, the analysis overpredicts AU^/AUrp (fig. 20). The average deviation 
between analytical and experimental values of AU^/AUj is 6. 13 percent for the diffuser 
type injectors compared to 17. 29 percent for the straight pipe injector. One of the 
assumptions made in the analytical program is that no radial temperature gradients exist 
in the ullage. This assumption is approximately true during expulsion with the diffuser 
type injectors (resulting in a fair analytical prediction of wall heating) but not true when 
using the straight pipe injector. Figure 21 compares the actual radial temperature pro- 
files and the predicted ones while using the straight pipe injector. 

Figure 22 presents a comparison of the AU^/AUrp ratio for different expulsion 
times for the three injector geometries. In all cases, between 12 and 20 percent of the 
total energy added to the tank appears as an increase in liquid energy. The data show no 
significant trend in AU^/AUp, with expulsion for any of the injector geometries. The 
probable errors associated with experimentally determining AU^ (see table II) are large 
enough to discourage any quantitative discussion, The straight pipe injector indicates 
less liquid heating (AU L , table II) than the diffuser type injectors. 

The decrease liquid heating and the increased evaporation obtained when using the 
straight pipe injector may be due to the inability of the liquid hydrogen to conduct heat 
from the surface to the bulk fluid at a rate equal to that at which the gas is adding heat to 
the surface. As a result, a local area of liquid at the surface quickly reaches the satur- 
ation point and evaporates, carrying with it some of the liquid's energy. The decrease 
in liquid heating when using the straight pipe injector may not occur for all tank geometries 
and sizes. 

The analytical data presented in figure 22 indicate an increase in liquid heating 
(increasing AUj/AU t ) for increasing expulsion times for all injectors. There is no 
appreciable variation in the predicted values of AU-p/AU™ for the three injectors at the 
various expulsion times . What little scatter there is is due to variations in initial conditions 
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Figure 21. - Comparison of experimental and analytic radial gas tem- 
peratures at various axial positions in the tank at end of 276-second 
expulsion. Straight pipe injector; inlet gas temperature, 289 K 
(520° R); tank wall thickness, 0. 762 centimeter (0. 30 in. ). 
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Figure 22. - Comparison of energy gained by liquid to total energy added ratio as function 
of expulsion time for three injector geometries. Inlet gas temperature, 284±6 K 
(511°±11° R); tank wall thickness, 0.762 centimeter (0.30 in.); tank pressure, 

34. 47x10^ newtons per square meter (50 psia). Curves faired through experimental data. 
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(e.g. , initial wall and gas temperature profiles, inlet gas temperature, and pressure) for 
the different injectors. 

The amount of usable propellant is an important parameter when considering the 
overall optimization of a pressurization system. For a pressure-fed system, all the 
propellant would be usable. However, for a pump-fed system, the amount of usable pro- 
pellant would depend on the temperature stratification buildup in the liquid and the net 
positive suctionhead (NPSH) of the pump. For a given outflow rate (expulsion time), the 
final temperature profile in the liquid was approximately the same for the various inlet 
gas temperatures and injector geometries used in this report. This means that for a 
given NPSH, the usable propellant would be the same regardless of inlet gas temperature 
or injector geometry. This result may not be true for different tank sizes or geometries 
or for different external heating rates. 

Figures 23 and 24 compare experimental and calculated ullage gas and wall temper- 
ature profiles for a diffuser type injector (radial) and the straight pipe injector at the 
end of approximately a 276 -second expulsion. 



Figure 23. - Comparison of analytical and experimental gas and wall temperatures at the end of 277-second ex- 
pulsion. Radial injector; inlet gas temperature, 288 K (518° R); tank wall thickness, 0. 762 centimeter 
(0.30 in.); run 486. 


38 



& 

£ 


ea 



Figure 24 - Comparison of analytical and experimental gas and wall temperatures at end of 276-second expul 
sion. Straight pipe injector; inlet gas temperature, 289 K (520° R); tank wall thickness, 0.762 centimeter 
(0.30 in.) ; run 492. 


The analysis overpredicts the tank wall temperatures (fig, 23) for the upper half of 
the tank wall when using the radial injector. The temperature over prediction accounts 
for the overprediction of 8. 03 percent on AU^y/AUrp (fig. 20). The analysis also over - 
predicts the ullage gas temperatures in the upper 70 percent of the ullage. In spite of 
this overprediction of ullage temperature, the analysis predicts the Mj/Mq ratio to 
within 1. 00 percent of that which was determined experimentally. This again is due to 
the ability of the analysis to predict the mass average temperature which was previously 
discussed. 

Figure 24 shows the effect that the straight pipe injector has on the wall and ullage 
gas temperature profiles. The experimental data show a nearly constant ullage temper- 
ature profile in the upper 40 percent of the ullage. The temperature then decreases al- 
most linearly to the liquid surface. The experimental wall temperature profile is similar 
to that of the ullage. The analytical gas temperatures increase almost linearly from the 
liquid surface to the top of the tank. There is good agreement between the analytical and 
experimental wall and gas temperature profiles in the lower 60 percent of the tank, but 
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the analysis overpredicts both gas and wall temperatures in the upper 40 percent of the 
tank. The predicted (constant radial temperatures) and experimental radial gas temper- 
ature profiles for the run data shown in figure 24 are compared in figure 21. The pre- 
dicted gas temperatures near the tank wall are with one exception higher than the exper- 
imental data. Thus, the larger overprediction of wall temperatures in the upper 40 per- 
cent of the tank coupled with the fact that the specific heat of aluminum rapidly increases 
with increasing temperature leads to the overprediction of 19. 42 percent on AlJ^/AlLj,. 
The good agreement between the analytic and experimental temperatures in the lower 
60 percent of the ullage (where approximately 87 percent of the total ullage mass is con- 
centrated) is the reason there is good agreement between the analytic and experimental 
mass ratios based on the net mass added experimentally (Mqj^). 

The experimental results, when using the straight pipe injector, indicate non-one- 
dimensional flow and substantial mass transfer, both of which are contrary to the 
analytical assumptions. The close agreement between analysis and experimental mass 
ratios (Mj/Mq N ) apparently is the result of one faulty analytic assumption compensating 
for another. These results would not necessarily be expected for tanks of different size 
or geometry or under different operating conditions. 

Moreover, the analysis does not predict the distribution of the energy losses (i.e. , 
AU^/AU^, or AU L /AU T ). It is generally concluded that the present analytical model is 
not strictly usable in its present form for the analysis of the straight pipe injector. 

Effect of tank wall thickness on pressurant requirements. - The results of refer- 
ence 3 indicated that between 77 and 93 percent of the total energy lost by the pressurant 
gas was lost to the tank wall. It was expected, therefore, that tank wall thickness would 
play an important roll in determining the pressurant requirements during the expulsion 
period. Additional tests were performed using a 1. 52-meter- (5-ft-) diameter spherical 
tank similar to the one already discussed except that the tank wall was chem -milled down 
to an average thickness of 0. 409 centimeter (0. 161 in. ). Although the chem -milling 
reduced the wall thickness by 46. 3 percent, the total weight of the tank was reduced by 
only 29. 5 percent due to the structural requirement of keeping the same thickness for the 
tank neck, lid, and girth support. 

Figure 25 presents a comparison of the Mj/Mq ratio for different expulsion times 
for the two tank wall thicknesses (0. 762 and 0. 409 cm; 0. 30 and 0. 161 in. ). Each tank 
used the hemisphere injector. The average inlet gas temperature for the 0. 409 centi- 
meter (0. 161 in. ) tank was 290 K (~520° R). The average inlet gas temperature for the 
0. 762 centimeter (0. 30 in. ) tank was 273 K (~490° R). The 0. 409 centimeter (0. 161 in. ) 
thick tank has an average of 6. 5 percent less pressurant requirement, higher Mj/Mq 
ratio, than the 0. 762 centimeter (0. 30 in. ) thick tank for all expulsion times. In most 
cases, the analysis again overpredicts the pressurant requirements (underpredicts the 
Mi/Mg ratio). See table in for deviations between analytical and experimental results. 
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.5— O .409 (.161) 1 to 5 

Open and closed symbols denote experi- 
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Total expulsion time, sec 

Figure 25. - Comparison of ideal pressurant requirement to actual pressurant require- 
ment ratio as function of expulsion time for two wall thicknesses. Hemisphere 
injector; inlet gas temperature, 281±15 K (506°±27° R); tank pressure, 34. 47X10 4 newtons per 
newtons per square meter (50 psia). Curves faired through experimental data. 


The 6. 5-percent decrease in pressurant requirements was less than expected. How- 
ever, the test with the 0. 409 centimeter (0. 161 in. ) thick tank had an average 17 K (31°R) 
higher inlet gas temperature. When based on the results already obtained on the effect of 
inlet gas temperature on pressurant requirements (fig. 8), this 17 K (31° R) higher inlet 
temperature would normally decrease the Mj/Mq ratio by approximately 4. 8 percent. 
Thus, for equal inlet gas temperatures, a nominal decrease in pressurant requirements 
of approximately 11. 0 percent would be expected for the 0. 409 centimeter (0. 161 in. ) thick 
tank. 

Figure 26 is a comparison of AU^/AUrp for different expulsion times for the two 
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Total expulsion time, sec 


Figure 26. - Comparison of energy gained by wall to total energy added ratio as function 
of expulsion time for two wall thicknesses. Hemisphere injector; inlet gas tempera- 
ture, 281±15 K (506°±27° R); tank pressure, 34.47X10 4 newtons per square meter 
(50 psia). Curves faired through experimental data. 
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tank wall thicknesses. The figure indicates that the same percentage of the total energy 
added to the tank was gained by the tank wall for both tank thicknesses. The absolute 
value of wall energy gained (table n) by the 0. 409 centimeter (0. 161 in. ) thick tank was 
approximately 6. 3 percent less than the 0. 762 centimeter (0. 30 in. ) thick tank. 

The absolute values of increase in ullage energy, liquid heating, and mass transfer 
are approximately the same (within measurement uncertainty) for the two tank wall 
thicknesses (see tables I and n). Thus, the 6. 5-percent reduction in pressurant require- 
ments when using the 0. 409 centimeter (0. 161 in.) thick wall tank is almost entirely due 
to the 6. 3 percent reduction in the energy gained by the tank wall. 

Figure 27 is a comparison of experimental final wall temperature profiles for the 
two tank wall thicknesses at the end of approximately 276-second expulsions. The 0. 409 
centimeter (0. 161 in. ) thick tank has much higher wall temperatures than the 0. 762 centi- 
meter (0. 30 in. ) thick tank in the upper position of the tank. The reduction in wall mass 
is nearly compensated for by the increase in wall temperature resulting in a difference 
in wall energy gain for the two profiles of only 6. 3 percent. The percent of the total wall 



50 100 150 200 

Wall temperature, K 


Figure 27. - Comparison of final wall temperature profiles at end of 
approximately a 276-second expulsion for two tank wall thicknesses. 
Hemisphere injector; inlet gas temperature, 278±7 K (500°±13° R); 
tank pressure, 34.47X10 4 newtons per square meter (50 psia). 
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energy gained by the tank lid for the two tanks was approximately the same (an average of 
21 percent of the total was gained by the lid for the thinner wall tank compared with 
24 percent gained for the 0.762 cm (0.30 in.) thick tank) even though the lid, for 
the thinner tank, represented a greater percent of the total tank weight (30 percent com- 
pared to 22 percent for the thicker tank). The reason for this was that the initial lid 
temperature for the thinner tank was an average of 30 K (54° R) higher than the thicker 
tank prior to expulsion. 

For the range of expulsion times investigated herein, the reduction of the wall thick- 
ness (tank mass) of 29. 5 percent resulted in lowering the pressurant gas requirement by 
6. 5 percent (11 percent when corrected for inlet gas temperature). For long duration 
expulsions where the final wall and ullage gas temperatures are dependent on inlet gas 
temperature and not tank wall mass, the reduction in pressurant requirements would be 
nearly proportional to the reduction in tank mass (due to a reduction in wall heat capacity). 

Effects of ramp rates, inlet gas temperature, and initial ullage v olume on the pres- 
sur ant require d to pressurize the tank from atmosphere to 34. 47xlQ 4 newt ons per square 
meter (50 psia) . - The amount of pressurant gas needed to initially pressurize a propel- 
lant tank may be important for certain missions. This is particularly so for multiburn 
missions where the tank is vented after each burn or where the coast period between 
firings is long enough to enable the ullage gas to collapse. 

As stated in the INTRODUCTION, the purpose of this investigation was to determine 
the capability of the analysis to predict the pressurant requirements during the initial 
pressurization period as well as the expulsion period. For this purpose, data were col- 
lected during the initial pressurization period for various pressurizing rates, inlet gas 
temperatures, and ullage volumes. 

Figure 28 is a comparison of the Mj/Mq ratio as a function of inlet gas temperature 
for two ramp rates at an initial ullage volume of approximately 4 percent. The data were 
taken using the hemisphere injector. At constant inlet gas temperature, the pressurant 
requirements decreased for increased ramp rates. The Mj/Mq ratio decreases for 
increasing inlet gas temperature. The reduced pressurant requirements (for the ramp 
period) for fast ramp rates using low inlet gas temperature is in agreement with the trends 
during the expulsion period. 

The modification of the analysis of reference 1 for the ramp period is discussed in 
appendix C . As can be seen in figure 28, the analysis is not capable of accurately pre - 
dieting the pressurant requirements during the initial pressurization of the 4-percent 
ullage. However, the prediction of the total pressurant requirement (initial pressuriza- 
tion and expulsion) is still good because the amount of gas required to initially pressurize 
the 4-percent ullage was only 4. 0 percent (maximum) of the pressurant requirement 
during expulsion. 
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Inlet gas temperature, °R 
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Figure 28. - Comparison of ideal pressurant requirement to actual pressurant requirement ratio as function 
of inlet gas temperature for two ramp rates at approximately 4-percent ullage volume. Hemisphere injector. 
Curves faired through experimental data. 
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Figure 29. - Comparison of ideal pressurant requirement to actual pressurant requirement ratio as 
function of ramp rate for various initial ullage volumes. Inlet gas temperature, 167 K (301° R). 
Curves fairedthrough experimental data. 
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The absolute perssurant requirements (both experimental and analytical) for the data 
presented in figures 28 and 29 are given in table IV. 

The comparison of Mj/Mq as a function of ramp rate for various initial ullage 
volumes at an inlet gas temperature of approximately 167 K (301° R) is shown in figure 29. 
This figure indicates increased Mj/M g for increased ramp rates for any given initial 
ullage volume. There is an increased Mj/Mq for the larger ullage volumes for a given 
ramp rate. Here again, the analysis is not capable of accurately predicting the pres- 
surant requirements. However, the analytical predictions seem to improve for the 


TABLE IV. - COMPARISON OF EXPERIMENTAL AND ANALYTICAL VALUES OF 
PRESSURANT GAS REQUIREMENTS FOR RAMP PERIOD 


Run 

Initial ullage 
volume, 
period 

Ramp rate 

Mass added 

Inlet gas 
temperature 

psi/sec 

(N/m 2 )/sec 

Experimental 

Analysis 

°R 

X 

ib 

kg 

lb 

kg 

14 

8.5 

2. 

54 

17.5X10 3 

0. 108 

0. 049 

0.088 

0.040 

524 

291 

17 

4.8 

2. 

54 

17.5 

.068 

.031 

.049 

.022 

526 

292 

60 

3.3 

1 . 

15 

7.93 

.068 

.031 

.048 

.022 

680 

378 

61 

3.4 

1 . 

15 

7.93 

.075 

.034 

.047 

.021 

720 

400 

64 

3.6 

2. 

54 

17.50 

.048 

.022 

.037 

.017 

723 

402 

65 

3.7 





.042 

.019 

.036 

.016 

726 

403 

67 

5.0 





.074 

.034 

.059 

.027 

715 

397 

68 

4.0 





. 057 

.026 

.036 

.016 

723 

402 

69 

1 

4.4 

1 . 

15 

7 

93 

. 086 

.039 

.057 

.026 

504 

280 

70 

4.0 

1 . 

15 

7 

93 

. 062 

.028 

.053 

.024 

336 

187 

71 

3.3 

2. 

54 

17 

50 

. 087 

.039 

.059 

.027 

335 

186 

73 

4.7 

2. 

54 

17.50 

.078 

. 035 

.066 

.030 

313 

174 

74 

3.0 

2. 

54 

17 

50 

. 073 

.033 

.056 

.025 

282 

157 

75 

3.2 

1 . 

15 

7.93 

. 086 

.039 

.049 

.022 

327 

182 

110 

28.2 

1 . 

48 

10.20 

.365 

.166 

.286 

. 130 

284 

158 

112 

27.8 

1 . 

48 

10.20 

.302 

.137 

.298 

. 135 

287 

159 

127 

55.4 

1 . 

15 

7 

93 

. 593 

.269 

.436 

. 198 

314 

174 

128 

54.8 

.66 

4. 55 

.844 

.383 

.520 

.236 

314 

174 

139 

55.6 

1 . 

48 

10.20 

.559 

.254 

.486 

.221 

307 

171 

380 

73.5 

1 . 

15 

7.93 

' .628 

i 

.285 

.550 

.250 

309 

172 

381 

75.4 

1 . 

48 

10.20 

i .622 

.282 

.601 

.273 

287 

159 
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larger ullage volumes. The transient process that occurs during the initial pressuriza- 
tion of the tank is too complex to be described by the present analytical model. The 
analytical program can, however, at least be used to predict the approximate magnitude 
of pressurant requirements during the ramp period. 

Concluding remarks . - The comparisons between the analytical and experimental 
results indicate that for the range of test conditions used, the analytical program and 
assumptions are adequate to allow prediction of gas requirements during expulsion when 
using a diffuser type injector. However, the analytical results are not as good when a 
straight pipe injector is used. The discrepancy is due to substantial amounts of evapor- 
ation and to the relatively large radial temperature and velocity gradients existing in the 
ullage during expulsion when a straight pipe injector is used. The analytical results are 
very good if compared to the net mass added to the ullage during the expulsion period for 
all injector geometries. 

The analysis is not capable of accurately predicting the pressurant requirements 
during the initial pressurizing period particularly for small initial ullage volumes. How- 
ever, since the absolute pressurant requirement for this ramp period for small initial 
ullage volumes is only a small fraction of the total pressurant required for expulsion, 
the analysis is considered good enough for approximation. 


SUMMARY OF RESULTS 

Tank pressurization and propellant expulsion tests were conducted to determine the 
effect of various physical parameters on the pressurant gas requirements. The experi- 
mental results were also compared with predicted results obtained from a previously 
developed analytical program . The analytical program was revised and extended for these 
tests conditions. The tests were conducted using two 1. 52-meter- (5 -ft-) diameter tanks 
with 0. 762 and 0. 409 centimeter (0. 30 and 0. 161 in. ) wall thicknesses. The tests were 
conducted at various inlet gas temperatures, outflow rates, and injector geometries with 
the following results . 

(1) The experimental results indicate decreasing pressurant requirements but 
decreasing Mj/M^ as the inlet temperature of the pressurizing gas is increased. Of all 
test variables, inlet temperature had the strongest effect on pressurant gas requirements. 

(2) The straight pipe injector had lower pressurant requirements and a greater 
M/Mg than the diffuser type injectors as a result of decreased wall heating, slightly 
less liquid heating, and increased liquid evaporation. Injector design had the second 
strongest influence on pressurant gas requirements. 

(3) Decreasing the tank mass by 29. 5 percent (by decreasing the tank wall thickness 
from 0. 762 to 0. 409 cm (0. 30 and 0. 161 in. )) decreased the pressurant requirements, 


46 



increased Mj/Mq, by an average of 6. 5 percent (11. 0 percent when corrected for varia- 
tion in inlet gas temperature). 

(4) The effects of inlet gas temperature and ramp rate on the pressurant gas required 
during the initial pressurization of the tank at a 4-percent ullage were as follows: 

(a) Increasing inlet gas temperature decreased the pressurant requirement and 
the Mj/Mq ratio for constant ramp rates. 

(b) Increasing ramp rate decreased the pressurant requirement and increases 
the Mj/Mq ratio for constant inlet gas temperatures. 

(5) The trends during the pressurization period for various initial tank ullage volumes 
and ramp rates are as follows: 

(a) Increased pressurant requirement but increased Mj/Mq for larger initial 
ullage volumes for a constant ramp rate. 

(b) Decreased pressurant requirements and increased Mj/Mq for increased 
ramp rates at a given initial ullage volume. 

(6) The ability of the analytical program to predict the experimental results was as 
follows: 

(a) The actual pressurant requirements were predicted within an average of 

12. 38 percent for all runs (maximum, 21. 28 percent; minimum, 6. 33 percent). How- 
ever, the analysis was able to predict the net mass added to the ullage to within an 
average of 4. 93 percent for all runs (maximum 16. 7 percent; minimum, 0 percent). 

(b) Tank wall heating was predicted within an average of 8. 04 percent for all runs 
with the exception of the straight pipe data which was within 17. 29 percent. 

(c) Liquid heating was predicted within a minimum of 12. 98 percent and a maxi- 
mum of 31.7 percent for all test conditions. 

(d) The analysis failed to predict accurately the pressurant requirements during 
the initial pressurization of the tank; however, it is considered to be good enough for 
approximation purposes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 7, 1969, 

180-31-02-01-22. 
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APPENDIX A 


VARIABLE GEOMETRY, HEAT LOSS TO TANK WALL, AND INTERNAL HARDWARE 

The basic analysis used in this report for predicting pressurant gas requirements 
was developed by W. H. Roudebush in reference 1 for a cylindrical tank. 

The major assumptions in the analysis of reference 1 are as follows: 

(1) The ullage gas is nonviscous. 

(2) The ullage gas velocity is parallel to the tank axis and does not vary radially or 
circumferentially . 

(3) The tank pressure does not vary spatially. 

(4) The ullage gas temperature does not vary radially or circumferentially. 

(5) The tank wall temperature does not vary radially or circumferentially. 

(6) There is no axial heat conduction in either the gas or the wall. 

(7) There is no mass transfer (condensation or evaporation). 

(8) There is no heat transfer from the pressurant gas to the liquid. 

Experiments performed at Lewis (ref. 3) confirmed most of these assumptions. The 

experimental results indicated, however, that there is significant heat transfer from the 
gas to the liquid with resulting mass transfer. 

For the purposes of this report, the analysis of reference 1 was modified for appli- 
cation to arbitrary symmetric tank shapes, and an attempt was made to incorporate the 
heat transfer from the gas to the liquid. The treatment of internal hardware (e.g. , tank 
baffles, instrumentation) was also modified to correspond to the treatment of heat transfer 
to the tank wall. 

The primary equations which deal with the pressurizing gas upon entering the tank 

are 

(1) The first law of thermodynamics 

(2) The continuity equation 

(3) The equation of heat transfer for a point in the tank wall 


First Law of Thermodynamics 


The form of the first law of thermodynamics used in the analysis in reference 1 for 
cylindrical tanks is 


9T 

at 


2h ZRT 

( T w 

rMPC w 
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- T) - V — + 

dx 


RTZj 


9P [ RTZq H C H 
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Modifying this equation to account for both arbitrary symmetric tank shapes and internal 
tank heat sinks gives 


9T 

at 


2h ZRT 

(T 

rMPC p w 



va T [ RTZ i ap | q h 
ax Mpc p at CpM H 


(Al) 


The first term on the right includes the effect of wall curvature. The last term, the 
energy lost to the internal hardware, is treated as the summarion of hardware compo- 
nents: (1) laminated thermoplastic, (2) stainless steel, and (3) copper. For the tanks 
in this investigation, 


_ A l h c< T l - T > , A 2 V T 2 - T > , A 3 h c< T 3 - T > 
M h pv pv pv 


(A2) 


where A j, Ag, and Ag are effective areas of materials contained within the volume 
element. 

Gluck and Kline, in reference 5, employed the free convection correlation to the 
pressurant gas (hydrogen, helium) for the pressurized transfer of liquid hydrogen: 

= Nu = 0. 13(GrPr) 1//3 (A3) 

k 


This correlation is used herein even though it was developed for cylindrical tanks . 
Pressurant gas transport properties were evaluated at the mean of the gas and wall tem- 
peratures. 


Continuity Equation (Area = f(x)) 

The basic form of the continuity equation for a cylindrical tank is presented in refer- 
ence 1 (eq. (24)) as 


9V _ Z 1 / 0T + V3T\ Z 2 9P 

ax zt \at 9x / zp at 

The modified form of the continuity equation due to variations in tank radius with 
distance along the vertical axis becomes 
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(A4) 


9V _ Z 1 /3T + V9T\ Z 2 3P 2V 0r 

ax zt \at ax / zp at r ax 

where Z ^ and Z 2 are defined in reference 1 as 



The last term in equation (A4) evolves from the derivation as follows. For the one- 
dimensional expression for continuity, 

— (pVA) + JL (pA) = 0 

ax at 

9 

The substitution A = irr is made where r is the position radius at location x along the 
vertical axis: 


-1- (pV r 2 ) + A (pr 2 ) = 0 

ax at 


The expression for density from the equation of state p = MP/ZRT is substituted: 



= 0 


The following velocity equation is obtained after performing the partial differentiation 
and after rearranging terms: 


av 

ax 



l ap _ 2v ar 
p at r ax 


When the expressions involving Zj and Z 2 are substituted in this equation, equation 
(A4) is obtained. 
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Tank Wall Heat Transfer 


Reference 1 (eq. (18)) gives the heat-transfer equation which represents the change 
in wall temperature as a result of the convective process for a cylindrical tank: 


at Ip C 
w^w w 

where q,„ is the rate of heat addition per 
w 

tank. 

For a small element of volume in the 


(T - T ) + — ^ (A5) 

1 w^w^w 

unit area to the tank wall from outside the 
x-direction, equation (A5) can be written as 


9T. 


W^r = h c A ^ T -V + Qw 


at 


(A6) 



For a wall of arbitrary shape, the following is evident from the sketch: 

A _ 2 nr As _ As _ _1_ 

V 2nra Ax a Ax Z w 

Therefore, equation (A5) holds also for this case. 

To account for the large mass concentration at the top of the tank, an equivalent l w 

was used. This Z was obtained by dividing the mass of the tank lid and flange con- 
w 

cections by the surface area of the first net point. 
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APPENDIX B 


EQUATIONS OF HEAT AND MASS TRANSFER AT THE GAS-LIQUID INTERFACE 

The energy and continuity equations (Al) and (A4) should be modified to incorporate 
both heat transfer from the ullage gas to the liquid surface and mass transfer into the 
analysis . 

The energy equation should incorporate two additional terms: 

(1) The heat-transfer rate from the ullage gas to the liquid interface (q^g) 

(2) The energy associated with mass transfer (M t A) 

Also, the boundary condition at the liquid surface should be revised to account for mass 
transfer. These additional terms can be related by performing an energy balance at the 
gas-liquid interface as done by W. A. Olsen in reference 6. The resulting relation, 
when the interface is saturated, is given by 


M t 

q U-S " q S— L + X ^ 

M 

= q S— L + X ^VAPOR" H L> ( B1 ) 

However, mass transfer was neglected since experimental data, presented in this report, 
indicate that the mass transfer rate was relatively small for most cases (and hence, 
the energy associated with the mass transferred is small in comparison to qg_X • The 
assumption was made that 


M. 


J.V. Lj j 

q U-S = q S— L + X ^VAPOR ■ h L* - ^ ( U l) 


(B2) 


The term d/dt(U L ) can be determined from experimental data. However, for the 
purpose of the analysis q^g must be related to the ullage gas variables. This is done 
by the relation 


q U-S - h c, l/ T 5 T saP 


(B3) 
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The flow process is free convection flow of pressurant gas down the tank wall and then 
radially inward across the liquid surface. 

The term T ga ^. was defined by the thermodynamic assumption of local equilibrium 
for a pure system with relatively small gradients as the saturation temperature corre- 
sponding to the tank pressure. 

With regard to h c reference 7 developed an equation from boundary layer theory 
for forced flow across a horizontal, semi-infinite, constant temperature flat plate given 
by 


Nu = 




(B4) 


The velocity of the gas across the liquid surface in terms of the gas velocity Vq 
down a vertical wall is given in reference 8 as 


V L = 0. 0975 V G 


(B5) 


where V G> obtained by solving the integrated energy and momentum equations at the wall 
boundary, is given by 


V r = 1. 185 -H 

G pz 


Gr 


1/2 


11/2 


|l + 0. 494 (Pr) 2//3 J 
Combining equations (B5) and (B6) for V L and substituting into (B4) give 


Nu = 


h c,L L 0.226(Pr) 1 / 3 (Gr) 1 / 4 


1 + 0. 494 (Pr) 2 / 3 ] 


1/4 


(B6) 


(B7) 


Equation (B7) shows that the conductance h T across the gas-liquid interface is simi- 

w j Xj 

lar in form to the empirical relation for free convection flow given in reference 9 as 


h_ t L 

Nu = c>lj = 0. 14(GrPr) n 
k 


(B8) 
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This equation with a value of n = 1/ 3, although somewhat arbitrary, is used in this 
investigation. 

At this point, some choice of T g , which is consistent with the definition of h c L 
and fits the data for q^ T _ g , (i.e., d/dt(U L )) must be made. 

In reference 10, which involved the pressurization of hydrogen with a low mixing 
diffuser and no liquid outflow, the adiabatic compression temperature given by 
T a( j = T Q (P/P 0 )fr"^/y was used as the choice for T g . This relation gave good agree- 
ment between analytical and experimental mass flux results. 

However, for the conditions described in reference 10, appreciable condensation 
occurred. For greater ullage gas mixing (due to diffuser characteristics as well as the 
liquid outflow process), T g would be expected to be a higher value than T ad . Refer- 
ence 10 indicated that as T g increases these is a tendency toward evaporation - that is, 
away from the condensation results that occurred when the adiabatic temperature was used. 

For the work described herein, T g was evaluated using the relation 
qy_g = h c - T ga P = d/dt (U L ) for a specific case. The value d/dt (Uj ) was deter- 
mined from an experimental run where the expulsion time was approximately 295 seconds. 
For this condition T g was determined to be 1. 5 times the adiabatic temperature, or 
approximately 53 K (95° R), for a tank pressure of 34. 47xl0 4 newtons per square meter 
(50 psia). This value of T g was used for all comparisons since most of the experimental 
ullage gas temperature profiles indicated a change in slope around 56 K (101° R). 

Using equation (B8), as well as the values for T g and T gat discussed above, gave 
the final form of the equation used to evaluate the heat transferred from the ullage gas 
to the liquid interface as follows: 


iu-S = r (0- M)(GrPr> 1/3 (53 K - T sat ) 
L 


In order to incorporate liquid heating to the analysis, the term 

4u-S A L or 14>(GrPr) 1 / 3 (53 K - T sat )A L 


X n V p 


LXV 
n p 


(B9) 


(B10) 


must be added to the right side of equation (Al). 
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APPENDIX C 


RAMP ANALYSIS 

An application of the work reported in reference 1 is the prediction of mass of pres- 
surant for the ramp and hold period. A separate computer program which determines 
mass of pressurant as well as the tank wall energy requirements during the ramp and hold 
period is described herein. The same equations which deal with the expulsion period, 
as outlined in appendix A, are still applicable in the ramp and hold analysis to predict 
mass of pressurant and wall energy requirements. This analysis computes the gas tem- 
peratures in the ullage at any time during the pressure rise from the gas energy equa- 
tion. The corresponding gas velocities are computed from the equation of continuity. 

The iterative method to be described shows how convergence is achieved in the solution 
of the gas* energy and continuity equations. 

The predicted mass of pressurant is based on an integration of the volume elements 
in the ullage at the end of the ramp and hold periods. 

Because of the complexity of the equations involved in the iteration and because of 
the small amounts of pressurant required for the ramp when compared with the expulsion, 
no mass transfer or energy to the liquid is included in this analysis. Mass transfer re- 
quirements accounted for approximately 10 percent of the experimental pressurant re- 
quirements for the ramp period. Quantitatively, the entire mass of pressurant require- 
ments for the ramp period was less than the expulsion period by a factor of 30 to 60 when 
the initial ullage was 5 percent of the tank volume. 


INPUT DATA REQUIREMENTS 

For the solution to proceed, a set of boundary and initial conditions are required. 
These conditions, which are the same for the expulsion as well as the pressurization, 
are as follows: 

(1) At time t = 0, the values of gas temperature T and wall temperature T W as 

functions of x, the position within the ullage 

(2) On the boundary x = 0, the value of inlet gas temperature T as a function of 

time 

(3) At the liquid surface, the value of gas temperature T, wall temperature T W’ 

and velocity V as functions of time (Although movement of the interface has 
been noted during the ramp pressurization period, no significant effect on the 
programmed output was noted with the value of V = 0 at the interface. ) 
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(4) Tank pressure P, outside heating rate q w , and inside hardware heating rate q H 

as functions of time (Like the other initial conditions, the pressure P as a 
function of time or ramp pressure curve is defined by a discrete set of points 
which approximate a smooth curve. In regions or pronounced curvature, more 
points are needed for accurate definition than for linear portions . ) 

(5) Constant value of heat transfer coefficient h , or a correlating equation from 

c 

which h c may be evaluated at each net point from values of T, T^, and P 

(6) Tank radius as a function of axial distance down from the top of the tank 

(7) Tank wall material properties: density p w and specific heat c w( T w) 

(8) Tank wall thickness (average membrane plus weld area thickness) as a function of 

axial distance down from the top of the tank 

(9) Pressurizing gas properties: molecular weight M, specific heat C p (T), and 

compressibility factor Z(P, T) 

(10) Initial ullage height, total time of run, the number of net points in the initial 

ullage space 

(11) The initial time step At used in following the pressure rise as well as establish- 

ing the points of computation 

(12) If the hold period is to be included in the analysis, then the time for the end of 

the ramp must be specified 


APPLICATION OF BASIC EQUATIONS 


Reference 1 makes the substitution of T . from the finite difference form of (A6) 

Wj 1 

into the finite difference form of the first law. Rearranging gives a quadratic in the gas 
temperature T! where the prime refers to a step forward in time and the quantities with- 
out the prime are evaluated at the previous time step: 


T! 2 


"‘( 1+v ‘ ^ ‘ “0 "Vi* 


% At 

V^vAv C w/ 


T! - anvJ-^T! , + T. 


Ax 


i-1 


= 0 


(Cl) 


The quantity marked with the asterisk may be evaluated either at the beginning or the 
end of the time interval. 

A difficulty can arise when evaluating the gas energy equation expressed by the 
previous quadratic. This occurs when the heat transfer takes place from the wall into 
the ullage gas. For this situation, the solution of the continuity equation provided nega- 
tive gas velocities which made it impossible for equation (Cl) to converge on the real 
roots. 
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At the start of the ramp (immediately after filling the tank), the initial wall temper- 
ature distribution in the ullage is higher than the gas temperature distribution. This is 
brought about since the wall surface above the liquid is exposed to the ambient temper- 
ature. But the ullage gas temperature near the liquid interface is close to the saturation 
temperature at one atmosphere. 

The technique used when T . > T. involved a direct substitution. The finite dif- 

W 2 1 1 

ference form of equation (A4) is 



(C2) 


Substituting this value of V- for V* in equation (Cl) results in the following cubic 
equation: 



+ a * “ T, - a* " ♦ <!,) - Oj’Ti — b, 

1 Ax \Z ) i At 1 1 Ax 11 1+1 1 1 X Ax 


T i 




(C3) 


This cubic equation is solved for the gas temperature T- . 


ANALYTICAL PROCEDURE 

The analytical procedure uses a variable time increment At in following the pres- 
sure rise. With this technique, the iteration was stable over a range of inlet conditions 
and the results were consistent with the recorded data. The iteration proceeds in the 
following manner. A flow diagram is shown in figure 30 for reference. 
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Figure 30. - Flow diagram showing temperature-velocity iteration in energy and continuity equations. 


The initial velocity distribution is determined at time t = 0 by substituting the initial 
values of dT/dt from equation (Al) into equation (A4). For all the ramp runs encountered 
in this investigation, an initial time increment of 1 second proved satisfactory. 


TEMPERATURE CALCULATIONS FROM TOP TO INTERFACE 

Since values of V have been obtained at each net point at time t = Tj = 0, attention 
is turned to equation (C3) which is cubic in T! . During the iteration, the cubic equa- 
tion (C3) is solved for the gas temperature T| starting at the point Ng in figure 31. 

When this equation is first solved for the ullage temperature distribution, a value for 
T| + j is not available. A substituted value of T i + 1 proved to be satisfactory as an initial 
guess to get convergence. All other quantities in equation (C3) are available from the 
initial conditions. The value for T! j, the temperature at Nj, is known as a boundary 
condition. 

The solution for T| at Ng follows, and this procedure continues to calculate gas 
temperatures until the boundary at the interface is reached. The values for the corre- 
sponding wall temperatures are calculated using the finite difference form of equation (A6). 


VELOCITY CALCULATIONS FROM INTERFACE TO TOP 

Although the ullage temperatures are computed starting at the top (fig. 31), the 
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velocity equation (C2) is used to calculate the ullage gas velocity starting with the point 
N near the interface. The velocity at the interface N , the boundary value, is zero 

y z 

with no expulsion. 

The ullage gas velocity is calculated from point to point until the top of the tank is 
reached. The new velocities are used in equation (C3) along with previous values of 
Tt + i and the temperature distribution is redetermined. This process is continued until 
convergence is achieved over the entire ullage. The time is then advanced to t 2 and a 
new set of velocities is determined. 

COMPLETING THE SOLUTION 

With the new velocities at time t 2 , we evaluate equation (C3) again starting at point 
N 2 and terminating at the interface . A value for T! +1 is always available from the 
previous iteration, although a substituted value of T^ + ^ is used as the first value. 

The new values for T! at all the points for time t 2 are used to recompute the ve- 
locity distribution. This new set of velocities is then compared with the previous set 
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and convergence is assumed if the deviation is less than half of 1 percent for every 
velocity in the time set. A time step is then taken to tg. 

If convergence is not achieved after 40 iterations, the time step is reduced and the 
iteration process is reinitiated. Generally the reduction in time step becomes necessary 
only when there is a severe change in the slope of the ramp curve particularly in the early 
stages of the pressure rise. 

For the new time tg, the temperature T! in equation (C3) is determined from its 
converged value using the iterative method. This procedure continues to evaluate the gas 
temperatures and velocity distribution in the ullage for each time step taken in following 
the rate of pressure in the tank. 



Figure 32. - Tank pressure as function of time 
during initial pressurization period for 
run 139. 
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The initial gas velocity distribution used in solving equation (C3) for each new time 
t is obtained from the previous time as follows: 




(C4) 


This iterative procedure can be used for a constant pressure representing the hold 
period. However, for initiating the ramp, an actual pressure rise must be used. A 
typical example is shown in figure 32 . 
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